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ABSTRACT 


The  structural  survivability  of  shipboard  mast/antenna  systems  subjected  to  underwater 
explosion  can  be  “designed  in,”  through  the  determination  of  the  structural  dynamics  of 
the  mast/antenna  system.  This  thesis  details  the  specialized  application  of  accurate  and 
efficient  analytic  methods  for  the  structural  dynamic  design  analysis  of  shipboard 
mast/antenna  systems.  Investigated  herein  are  a  class  of  substructuring  methods,  generally 
referred  to  as  component  mode  synthesis  methods,  which  provide  for  the  rapid  calculation 
of  dynamic  response  of  the  mast/antenna  structural  system  to  weapons  effects. 
Additionally,  the  methods  also  provide  for  the  simulation  of  live  fire  testing.  The  methods 
allow  the  individual  antennae  and  the  mast  ea^  to  be  independently  modded,  arbitrarily 
combined,  and  the  combined  system  dyiuimic  response  rapidly  caloilated  to  determine  the 
structural  survivability  of  a  proposed  mast/antenna  configuration.  This  rapid  and 
‘‘modular^’  component-based  analysis  capability  is  specifically  tailored  for  interactive 
computer-aided  design  analysis  of  shipboard  mast/antenna  systems. 
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L  INTRODUCTION 


This  thesis  documents  the  development  of  analytic  methods  for  maximizing  the  combat 
survivability  of  shipboard  structural  systems  subjected  to  weapons  effects.  Survivability 
will  be  improved  through  the  characterization  of  the  mast/antenna  system  structural 
dynamics  and  the  development  of  specialized  deagn  analysis  tools  for  the  prediction  and 
minimization  of  dynamic  response  due  to  weapons  effects.  The  objective  is  improved 
system  combat  survivability. 

Additionally,  the  methods  will  be  developed  in  the  context  of  the  analytical  simulation 
of  live  fire  test  and  evaluation  (LFT&E)  for  shipboard  systems.  Those  shipboard  structural 
systems  which  undergo  linear  elastic  dynamic  response  due  to  live  fire  effects  can  be 
evaluated  for  live  fire  survivability  using  the  simulation  methods  to  be  developed,  therdjy 
eliminating  the  need  for  actual  LFT&E  for  these  systems.  Alternatively,  these  simulations 
will  be  of  benefit  in  the  planning  of  actual  live  fire  test  and  evaluation  (LFT&E)  programs. 
The  results  herein  focus  on  shipboard  mast/antenna  structures.  Shipboard  mast/antenna 
systems  must  be  designed  to  withstand  moderate  to  severe  shock  loading  induced  by 
underwater  explosion  (UNDEX)  of  conventional  or  nuclear  type.  The  UNDEX  delivers 
devastating  forces  to  the  targets  in  the  form  of  incident  shock  wave  pressure,  gas  bubble 
oscillation,  cavitation  closure  pulses,  and  various  reflection  wave  effects.  These  shock- 
induced  forces  then  propt^te  through  the  ship  to  the  various  systems,  equipment,  and 
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top-side  structures  including  the  mast  and  antennae.  The  response  of  the  mast  and 
antoinae  to  the  UNDEX  shock  wave  is  basically  linear  elastic  and  vibrational  ir.  nature. 
The  mast  and  antennae  tend  to  vibrate  at  their  fundamental  natural  frequency,  or  at  a  low 
range  of  natural  frequencies.  The  maximum  amplitude  of  the  vibration  usually  occurs  after 
the  shock  wave  passes  the  ship.  The  shock  response  wave  form  is  remarkably  different  at 
various  levels  within  the  ship.  In  essem^,  the  ship  acts  as  a  low  pass  structural  filter  which 
alters  the  characteristics  of  the  propagating  shock  wave  from  one  possessing  high 
frequency  components  to  one  that  contains  relatively  low  fiequency  components  [Ref  l;p. 
2].  Thus,  the  shock  survivabUity  of  the  mast/antenna  system,  which  is  located  top-side,  is  a 
vibration  problem  in  which  relatively  low  fiequency  equipment  support  excitations  are 
observed.  The  emphaas  on  design  analy^  relates  directly  to  the  survival  of  the  mission 
critical  systems  on  the  platform.  The  ability  of  the  naval  vessel  to  cany  out  its  misaon 
after  being  subjected  to  an  UNDEX  threat  depends  on  the  survivability  of  these  systems, 
and  specifically  the  mast/antenna  system.  Combat  survivability  of  new  systems,  such  as  the 
mast/antenna  system  can  be  “deagned  in”  by  accounting  for  the  structural  dynamics  of  the 
system  during  the  design  process.  The  methods  developed  herein  focus  on  the  the 
structural  dynamics  of  the  mast/antenna  systems,  so  that  their  combat  survivability,can  be 
directly  addressed  in  the  deagn  process.  Additionally,  the  methods  will  make  possible  the 
improvement  of  survivability  of  existing  systems.  For  example,  survivability  can  be 
improved  by  dynanucaUy  tuning  and  relocating  antennae  base^  on  the  application  of  the 
methods  to  be  described. 
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A.  BACKGROUND 


The  dynamic  response  of  a  shipboard  antenna  is  dependent  on  the  dynamic  interaction 
of  the  antenna  with  the  mast  during  response  to  weapons  effect.  Large  dynamic  loads  in 
an  antenna  can  result  if  (a)  the  antenna  is  mounted  on  the  mast  at  a  location  with  large 
accelerations  due  to  weapons  effects,  or  (b)  the  antenna  has  its  natural  frequencies  in  close 
proximity  to  the  excited  natural  frequencies  of  the  mast.  In  recent  years,  the  Navy  has  had 
frequent  occurrences  of  shipboard  antennae  systems  failing  structurally  after  being 
subjected  to  shock  due  to  weapons  effects  [Ref.  2].  In  order  to  design  these  structural 
systems  (i.e.  mast  and  antennae)  for  minimum  dynamic  response  and  hence  maximum 
survivability,  the  structural  dynamic  parameters  which  determine  the  dynamic  response  of 
the  system  must  be  accurately  quantified.  Ttw  primary  structural  dynamic  parameters  to  be 
determined  are  the  modal  parameters  (i.e.  natural  frequencies,  mode  shapes,  modal  mass 
and  damping)  of  the  mast  and  the  various  individual  antenna.  The  modal  parameters  are 
required  to  characterize  the  structural  dynamics  of  each  substructure,  e.g.  the  mast  and 
each  antenna,  and  hence  characterize  the  dynamics  of  the  combined  structural  system. 
Given  an  accurate  coupled  system  analytic  dynamics  model,  weiqrons-induced  dynamic 
re^nse  can  then  be  predicted,  and  system  designs  can  be  evaluated  and  optimized  with 
reflect  to  survivability.  The  coupled  system  analytic  dynamics  model  can  serve  as  the 
basis  for  the  computer  simulation  of  LFT&E. 
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B.  OVERVIEW  OF  THE  SUBSTRUCTURE  APPROACH  TO  THE  DESIGN 

ANALYSIS  OF  MAST/ANTENNA  SYSTEMS 

The  methods  described  herein  are  directed  at  the  automated  design  analysis  of 
mast/antenna  systems.  The  methods  provide  accurate  estimates  of  the  modal  parameters 
for  a  mast/antenna  structural  system,  and  therefore  will  provide  accurate  estimates  of  the 
dynamic  response  due  to  weapons  effects.  Generally  referred  to  as  **component  mode 
synthesis,”  these  substructuring  methods  make  use  of  independent  finite  dement  modds 
for  the  mast  and  each  antenna.  In  ord»  to  allow  a  designer  to  rapidly  assess  for 
survivability  a  large  number  of  candidate  mast/antenna  system  designs,  the  methods  are 
computationally  effident  as  wdl  as  accurate.  With  respect  to  mast/antenna  systems,  the 
component  mode  synthesis  process  will  allow  a  designer  to  analytically  “install”  the 
various  antenna  models  into  the  mast  model,  and  rapidly  calculate  coupled  mast/antenna 
S3rstem  UNDEX  dynamic  response.  Wl^n  incorporated  into  a  computo’-aided  design 
environment,  the  complexities  of  the  calculation  wUl  be  transparent  to  the  designer,  and 
will  allow  the  incorporation  of  sdf-checks  and  protection  against  user  error  and  misuse. 
The  substructure  approach  to  mast/antoma  structural  dynamic  analysis  can  be  briefly 
outlined  as  follows: 

•  A  designer  either  finds  the  dynamic  characteristics  of  the  various  antennae  to  be 
in^alled  from  a  “catalog”  (database)  of  antenna  modal  parameters,  or  calculates 
individual  antoma  modal  parameto^s  from  a  finite  element  model  of  the  antenna.  The 
modal  parameters  of  the  antenna  constitute  the  antenna  dynamic  model. 
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•  The  various  antennae  dynamic  models  are  analytically  coupled  with  the  mast  model, 
and  the  dynamic  re^nse  of  the  coupled  mast/antenna  sys«n  due  to  weapons  Elects 
is  calculated.  If  unacceptable  dynamic  response  levels  are  calculated,  the  various 
antennae  models  can  be  rapidly  repositioned  on  the  mast,  or  exchanged  with  other 
antennae,  and  the  new  dynamic  response  calculated. 

This  scheme  has  several  significant  advantages  for  the  automated  design  anal3rris  of 
mast/antenna  systems.  The  primary  advantages  include: 

•  The  ability  of  these  methods  to  treat  the  mast  and  antennae  as  “substructures,”  and 
arbitrarily  and  repeatedly  combine  them  for  the  rapid  calculation  of  dynamic 
response  will  make  possible  the  evaluation  of  a  greater  number  of  mast/antenna 
configurations,  and  hence  will  greatly  &cilitate  the  determination  of  an  optimal 
configuration  with  respect  to  combat  survivability. 

•  The  various  masts  and  antennae  are  fabricated  by  various  independent  contractors. 
The  conqranent  mode  synthesis  method  aUows  the  separate  modeling  of  the  mast 
and  antennae,  and  tho’efore  naturally  preserves  the  indq)endence  of  the  contractors. 

•  The  formulations  to  be  described  are  modal,  and  therefore  can  fimction  equally  well 
with  analytically  derived  modal  parameters,  or  with  modal  parameters  idoitified  in  a 
vibration  test. 

The  analytic  methods  for  the  genotudon  of  the  coupled  mast/antenna  model  are  the 
fiscus  of  this  work.  To  be  evaluated  in  this  report  are  sevoal  component  mode  methods 
for  substructure  synthesis:  the  Craig-Bampton  method  and  two  residual  flexibility 
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fonnulations.  The  methods  are  specialized  for  the  mast/antenna  analysis  problon,  and  their 
rdative  merits  ctunpared  in  the  context  of  combat  survivability.  The  methods  are !  ased  on 
the  modal  represoitation  of  components;  that  is,  rather  than  representing  a  structure  using 
the  mass  and  stiffiiess  matrices  generated  in  a  finite  element  model,  these  methods  employ 
various  classes  of  **mode  shapes”  to  represent  the  substructures  or  components.  For 
example,  the  &mitiar  normal  modes  of  vibration  are  one  class  of  mode  shape  used. 

The  computational  efficioicy  of  these  methods,  which  is  critical  to  their  efifectivraess 
in  a  con^Hiter  aided  design  environment,  comes  fi'om  their  ability  to  accurately  describe  a 
component  with  a  minimum  number  of  mode  shapes.  The  sections  of  this  report  which 
follow  will  describe  the  above  mentioned  synthesis  formulations,  aixl  demonstrate  their 
rdative  accuracy  and  e£5dency  in  the  calculation  of  the  dynamic  response  of  a  small  yet 
representafive  mast/antenna  modd,  subjected  to  a  variety  of  applied  harmonic  forces  as 
well  as  deck  accderations  and  displacements.  The  modd  used,  which  includes  a  mast  and 
a  single  antenna,  is  of  a  smaU  size  compared  with  that  required  to  represent  an  actual 
mast/antenna  structures.  However,  the  modd  has  all  the  features  necessary  to  allow  the 
assessment  and  critical  analysis  of  the  component  mode  synthesis  methods. 

Specifically,  the  three  synthesis  methods  will  each  be  used  in  the  following  analyses; 

(1)  Caknlation  of  ouut/anteniui  coa|ried  system  modal  parameters:  This  is  the 
fimdamental  assessment  of  a  method’s  accuracy.  Prior  to  performing  the  synthesis,  modal 
parameters  are  calculated  for  the  antenna  modd  and  the  mast  model.  The  appropriate 
component  representation  is  generated  and  the  niast/amenna  system  is  synthedzed.  The 
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coupled  system  natural  frequencies  are  calculated  and  are  compared  with  the  natural 
frequoicies  calculated  using  a  standard  finite  element  procedure.  The  standard  finite 
element  procedure  means  the  assembly  of  a  single  model  represeming  the  total 
mast/antenna  system.  A  comparison  of  floating  poim  operations  (FLOPS)  accumulated  in 
all  cases  is  also  provided.  This  comparison  will  demonstrate  the  computational  advantage 
of  the  synthesis  methods,  an  advantage  critical  to  the  developent  of  an  automated  design 
analysis  system. 

Using  the  synthesized  mast/antenna  modd,  the  following  analyses  are  presented; 

(2)  Calculation  of  antenna  peak  dis|rfacment  due  to  harmonic  forcing:  A  simple 
harmonic  forcing  function  is  ^plied  to  the  mast  and  the  peak  displacement  of  the  antenna 
free  Old  C*f>P’0  is  calculated,  again  using  all  three  component  mode  syntheds  methods,  as 
well  as  udng  a  standard  finite  element  procedure. 

(3)  Calculation  of  mast/antenna  interface  internal  stresses  due  to  harmonic 
forcing:  A  simple  harmonic  forcing  function  is  applied  to  the  mast  and  the  bending 
moment  and  shear  loads  in  the  mast/antenna  connection  are  calculated.  Note  that  these 
imemal  loads  are  directly  proportional  to  stress,  and  hence  are  the  critical  quantities  which 
must  be  calculated  in  order  to  assess  structural  survivability.  These  calculations  are 
repeated  for  all  three  sythesis  formulations,  as  well  as  for  the  standard  finite  element 
procedure. 
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IL  FORMULATION  OF  FINITE  ELEMENT  MODEL  AND  GENERAL 
COMPONENT  COUPLING  PROCEDURES 

The  theory  presented  herein  is  taken  directly  and  exclusively  from  reference  [12], 

As  discussed  in  the  Introduction,  the  finite  dement  (FE)  procedure  will  be  employed  to 
generate  mathematical  models  of  the  components  (substructures)  involved,  namdy  the 
mast  and  the  antenna.  The  FE  procedure  produces  stiffiiess,  mass,  and  less  commonly, 
damping  matrices  which  represent  the  structural  dynamics  of  each  components.  In  orda*  to 
fiuthfiilly  capture  the  geometric  and  material  conqjlexities  of  these  components,  the  finite 
element  discretization  must  necessarily  involve  many  degrees-of-fieedom  (DOF),  and 
hence  \ie  above  mentioned  system  matrices  can  be  quite  large.  The  time  and  cost 
associated  with  the  extraction  of  the  modal  paramters  (natural  frequencies,  mode  sh^)es, 
and  oKxlai  mass)  fix)m  these  large  matrices  precludes  the  performance  of  the  repeated 
design  analyses  required  to  arrive  at  an  optimal  design.  The  conq)onent  mode  synthesis 
methods  bypass  the  repeated  extraction  of  the  modal  parametm  for  a  complete 
mast/antenna  system  by  directly  using  the  modal  parameters  calculated  for  each 
component.  The  calculation  for  the  componoit  **modes”  is  performed  once  for  each 
component,  and  the  total  system  dynamics  are  synthesized  uang  the  various  sets  of  modes 
so  calculated.  The  ^nthesis  methods  not  only  provide  very  accurate  predictions  of 
^namic  response,  but  also  provide  a  substantial  decrease  in  the  time  required  to  compute 
dynamic  response,  hence  allowing  tl»  performance  of  additional  design  analysis  iterations. 
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A.  FINITE  ELEMENT  FORMULATION 

Although  FE  modeling  typically  involves  the  full  range  of  element  types  available  (e.g. 
beam,  plate,  shell),  for  purposes  of  this  thesis  the  antenna  and  the  mast  will  each  be 
modeled  using  beam  elements  only.  This  model,  although  simple,  is  all  that  is  necessary  to 
investigate  the  various  component  coupling  procedures.  All  methods  presented  herdn  are 
applicable  to  any  structural  model,  and  the  results  and  conclusions  presetted  are  directly 
applicable  to  the  analysis  of  structural  systms  of  any  complexity. 

Traditionally,  the  mast  and  antennae  are  modeled  together  as  a  system.  Alternatively, 
the  mast  and  antennae  can  be  modeled  separately.  By  modeling  the  mast  and  antennae 
separately,  several  benefits  arise: 

•  Masts  and  antennae  are  generally  fabricated  by  different  defense  contractors. 
Therefore,  modeling  the  mast  and  antennae  separately  would  best  preserve  this 
independence. 

•  Modeling  the  mast  and  antennae  separately  would  permit  the  development  of  a 
single  data  file  containing  only  mast  design  specifications,  and  several  separate  data 
files  containing  antennae  design  specifications,  one  datafile  for  each  amennae. 
With  this  modular,  component-based  ^proach  comes  the  flexibility  of  exchanging 
antennae  and/or  changing  antennae  placement.  This  allows  the  rapid  assessment  of 
many  mast  and  antennae  configurations  for  dynamic  response  characteristics. 

By  modeling  the  mast  and  antennae  separately,  the  computational  efficiency  increases  as 
compared  to  modeling  the  mast  and  antennae  together.  This  computational  advantage  is 
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due  to  tte  &ct  that  the  cost  associated  with  the  calculation  of  the  modal  parameters  for  a 
single  structural  modd  is  proportional  to  the  cube  of  the  number  of  DOF  of  th '  model 
[Ref  3:  p.  231].  The  calculations  performed  herein  demonstrate  this  comparison  between 
a  total  mast/antoma  model  and  a  modd  drived  from  the  synthesis  of  mast  and  antenna 
substructure  modds.  The  benefit  is  associated  not  just  with  the  calculation  of  the  modal 
parameters,  but  also  with  the  calculation  of  dynamic  response  to  assess  UNDEX 
survivability. 

B.  GENERAL  COMPONENT  COUPLING  PROCEDURES 

The  term  **component  mode  synthesis’^  refers  to  the  manner  in  which  each  substructure 
is  mathenuitically  represented  prior  to  coupling,  and  is  based  on  a  truncated  modal 
esqiansion.  This  representation  is  most  femiliar  in  the  context  of  the  calculation  of  dynamic 
response.  Here,  the  dynamic  re^nse  of  a  structure  can  be  written  as  a  linear  combination 
of  the  mode  sluq)es  calculated  for  the  structure.  If  the  fiequency  range  of  excitation  is 
contained  in  the  frequency  range  of  the  calculated  modes,  then  the  dynamic  response 
calculated  using  the  modes  will  be  of  acceptable  accuracy.  Of  course,  the  question  of  how 
many  modes  to  retain  is  non-trivial  and  prd>lem  specific.  However,  the  coiiq)utational 
efifidet^  of  a  modal  approach  to  structural  dynamics  including  the  component  mode 
synthesis  methods  to  be  presented,  comes  from  the  retention  of  a  number  of  modes  which 
constitute  a  mathonatical  model  much  smaller  than  the  original  mass  and  stififiiess  matrices 
firom  which  the  modes  were  calculated. 
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Component  mode  synthesis  makes  use  of  several  types  of  vibrational  mode  shapes, 
distinguished  by  the  boundary  conditions  imposed  on  the  substructure  prior  to  the 
calculation  of  these  mode  shapes.  In  addition  to  these  vibration  mode  shapes,  the  various 
component  mode  synthesis  methods  require  additional  types  of  mode  shapes  to  be 
calculated  and  included  with  the  vibrational  modes.  Therefore,  the  term  "component  mode 
synthesis”  (CMS)  is  a  suitable  name;  a  single  structure  is  synthesized  from  separate 
substructures  and  each  substructure  is  mathematically  represented  by  an  appropriate  set  of 
mode  shapes,  calculated  from  the  finite  elemoit  model  of  each  substructure.  The  following 
are  definitions  of  the  various  types  of  mode  shapes  that  are  used  in  the  component  mode 
synthesis  formulations  investigated  herein. 

1.  Free  Interface  Normal  Modes 

The  free  inter&ce  normal  modes  are  the  modes  of  the  component  when 
unrestrained  at  all  inters  DOF.  From  the  LIST  OF  SYMBOLS  AND 
ABBREVIATIONS,  the  interface  coordinates  are  denoted  by  the  subscript  "I"  and  the 
internal  coordinates  are  denoted  by  the  subscript  "O".  The  interface  coordinates  are  the 
coordinates  where  the  substructures  are  coupled.  The  internal  coordinates  are  all 
coordinates  that  are  not  interface  coordinates.  Free  interface  normal  modes  are  calculated 
by  solving  the  following  eigenvalue  problem; 


(1) 
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The  ^iffiiess  and  mass  matrices  in  Eq.  (1)  are  partitioned  as  follows: 


kJ  mJ 


The  number  of  equations  defined  by  Eq.  (1)  is  equal  to  the  number  of  rows  or 
columns  in  the  mass  and  stifihess  matrices.  The  number  of  columns  or  rows  in  [K]  or  [M] 
equals  the  number  of  DOF  of  the  component  in  physical  coordinates. 

2,  Fixed  Interface  Normal  Modes 

The  fixed  interfiice  normal  modes  are  the  modes  of  the  component  restrained  at  its 
interfece  DOF.  The  fixed  inter&ce  noimal  modes  have  the  following  form; 


The  upper  partition  of  the  fixed  inter&ce  normal  modes,  or  },  is  obtained  from  the 

solution  to  the  following  eigenvalue  problem; 


(2) 


In  words,  the  matrix  of  fixed  interfiu^  normal  modes  is  a  partitioned  matrix 
consisting  of  the  matrix  of  mode  shapes  obtained  in  the  solution  to  Eq.  (2)  in  the  upper 
partition,  and  a  matrix  of  zeros  in  the  lower  partition.  The  zeros  imply  zero  displacement 
at  the  inter&ce.  The  number  of  rows  in  the  matrix  of  zeros  is  equal  to  the  number  of 
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interfilce  coordinates,  while  the  number  of  columns  is  equal  to  the  number  of  internal 
coordinates  of  the  substructure. 

Both  the  upper  partition  of  the  fixed  inter&ce  normal  modes  and  free  interface 
normal  modes  are  unity  modal  mass  normalized  such  that  the  following  property  is 
satisfied; 


(3) 

[♦''r[Ml[<^N]=[I]  (4) 


3.  static  Constraint  Modes 

Static  constraint  modes  are  calculated  by  enforcing  a  unit  deflection  on  each 
interface  DOF  while  holding  all  other  DOF  restrained.  Calculating  the  resulting 
displacements  of  the  internal  coordinates  defines  the  static  constraint  modes.  Thus,  the  set 
of  static  constraint  modes  is  defined  bv  the  equation; 


Koo 

Ko.‘ 

^01  ■ 

po,l 

LK» 

bn  J 

LRoJ 

(5) 


where  [Rn]  is  the  numix  of  “reactions”  at  the  interface  or ‘T’  coordinates. 
From  the  top  row  partition; 


E^oi]  [*^00^01] 


(6) 
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If  h  can  be  assumed  that  no  external  forces  or  inertial  forces  are  applied  to  the  internal 
DOF,  as  in  a  static's  problem,  the  matiix  of  static  constraint  mode  sh^s,  or  [$^],  ar  3 
represented  by  the  following  equation; 


where  [<^^]is  the  static  constraint  mode  matrix. 

4.  Rigid  Body  Modes 

Rigid  body  modes  are  possessed  by  systems  that  are  not  restrained.  Rigid  body 
modes  have  zero  frequency.  They  can  be  solved  for  using  Eq.  (2).  They  can  also  be 
solved  for  in  the  same  way  that  static  constraint  modes  are  solved,  provided  that  the 
number  of  coordinates  retained  is  equal  to  the  number  of  rigid  body  modes.  The 
eigensolver  in  MATLAB™,  the  software  that  was  used  in  the  examples  in  Chapters  IV 
and  V,  is  ineflfective  in  producing  rigid  body  modes  using  Eq.  (2).  The  rigid  body  modes 
and  associated  fiequencies  as  calculated  by  MATLAB  using  Eq.  (2)  possess  a  complex 
part.  Therefore,  rigid  body  modes  can  be  obtained  using  Eq.  (S).  However,  it  was 
determined  when  applying  the  property  in  Eq.  (4),  the  orthogonality  property,  that  the 
rigid  body  modes  are  linearly  independent  but  not  orthogonal  with  respect  to  the  mass 
matrix.  To  produce  orthogonal  mode  shapes  from  a  set  of  mode  shapes  that  are  linearly 
independent  a  theory  fiom  linear  alg^ra,  the  Gram-Schmidt  theory  [Ref  4;  p.  I6S],  was 
invoked.  The  Gram-Schmidt  theory  will  now  be  presemed. 
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Given  three  linearly  independent  vectors  vectors  are 

selected  because  there  are  three  rigid  body  modes.  In  general,  six  vectors  would  need  to 
be  made  orthogonal  corresponding  to  six  ri^d  body  modes)  that  three  linearly  independent 
and  orthogonal  vectors,  ^  {^^2  {^3 }  obtained  as  follows; 

Find  a  vector  ^2 }  the  following  orthogonality  property  is  satisfied 

(8) 

To  satsify  the  requirements  of  Eq.  (8),  ^2 }  be  defined  by: 

where  a  is  a  scalar  that  is  used  to  extract  the  components  of  }  that  lie  in  the  vector 
space  of  {vj }. 

Substituting  Eq.  (9)  into  Eq.  (8)  the  following  relation  is  obtained; 

{v2  ^  (10) 

Expanding  Eq.  (10),  the  following  relation  is  obtained; 

{vj^  [M]^,}-a  [M]^,}=0  (11) 
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From  Eq.  (1 1)  a  is  obtained  as  follows: 


(12) 


By  substituting  Eq.  (12)  into  Eq.  (9),  the  second  linearly  independent  and  orthogonal 
vector^, }  ^  obtained. 

Now  there  are  two  linearly  independent  and  orthogonal  vectors  with  respect  to  the 
mass  matrix.  A  third  linearly  independent  and  orthogonal  vector  will  now  be  determined. 

Using  the  same  formulation  as  above,  a  vector  ]is  defined  as  follows: 

(13) 

where  0  and  y  are  scalars  that  are  used  to  extract  the  components  of  and 
respectively  that  lie  in  the  vector  space  of  {/, }. 

Uang  the  properties  of  orthogonality  the  following  relations  must  be  satisfied: 

(14) 

=  0  (15) 
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Substitutiiig  Eq.(13)  into  both  Eq.(14)  and  Eq.(lS)  the  following  relations  are  obtained: 


{v3  -/3vj  y  [M]^, }  =0  (16) 

^3  -/3vj  y  [M]^j }  =0  (17) 

By  expanding  Eq.(lb)  and  Eq.(17).  the  following  relations  are  obtained: 


Solving  Eq.(18)  and  Eq.(19)  simultaneously  the  following  expressions  are  obtained  forjS 
and  y: 

S  tMl^3}fr3y tM3^,}~{v,y [M]f 3} 

^2  f  [Ml{v,  }^,  7  [M]^,  }  -^3  J  [M]^3  }^,  f  [Ml^, } 

-^2  Y  [m]{v3  }f  3  y  [mi^,}+^3  y  [m]^,}{v3  y  } 

Substituting  Eq.  (20)  and  Eq.  (21)  into  Eq.  (13),  the  third  linear  independent  and 
orthogonal  vector  {9^3]  is  obtained. 

Therefore,  rigid  body  modes  will  be  solved  for  using  Eq.  (S)  and  the  Gram-Schmidt 
procedure  derived  above. 


(20) 


(21) 
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5.  Rctkiaal  Ftaubflity  Modes 

Before  defining  "residual"  flexibility,  the  concept  of  flexibility  must  first  be  «!efined. 
The  fledbility  of  a  restrained  structure  (i.e.  a  structure  whose  stifihess  matiix  is  of  full 
rank)  is  the  inva%  of  the  stiffiiess  of  the  structure.  By  inverting  the  stiffness  matrix,  one 
obtains  the  flexibility  matrix  as  follows; 


[G]  =[K]‘'  (22) 

Equation  (22)  can  also  be  wiittco  as  fi>llows: 

[Gl  -fb"tA  J’^"r  (23) 


The  residual  flexibility  matrix  is  obtaimd  from  the  flexibility  matrix,  the  kept  free  interfiice 
normal  modes,  and  the  inverse  of  the  natural  frequencies  as  follows: 


[o']  -[G]  K*D"rAD,]"[*D”r  («) 


The  residual  flexibility  modes  are  the  portion  of  the  exact  static  flexibility  shapes  that  are 
not  rqjresented  by  a  set  of  retained  modes.  Residual  flexibility  modes  require  the 
knowledge  of  other  modes  that  are  retained  in  the  model  and  are  dependem  upon  the 


retained  modes.  There  are  two  ways  to  calculate  the  residual  flexibility  modes. 

1)  If  the  structure  is  grounded,  such  as  the  mast,  then  the  stiffiiess  matrix  is  fliU 


rank  and  invertible.  By  post-multiplying  Eq.  (24)  by 
flexibility  modes  for  restrained  substructures. 


Jn . 


one  obtains  the  reudual 
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W  “{O'']- 

Lrn  . 


(25) 


2)  However,  If  the  component  is  not  grounded  before  assembly,  such  as  an 
amenna,  then  an  inertia  relief  solution  must  be  calculated  to  determine  the  flexibility  matrix 
as  follows; 


[G]  (26) 

where  i'*  are  the  rigid  body  nuxles  of  the  structure. 

[K  *]is  formed  by  inverting  the  restrained  or  intanal  partition  of  the  stiffness  matrix  in  the 
following  way; 


(K*l  = 


0  0 
?  Koo*’. 


(27) 


This  "new"  [G]  or  flexibility  matrix  is  flee  of  rigid  body  modes.  The  Craig-Chang 
fonmilation,  which  will  be  presented  in  Chapter  III  uses  firee-interface  normal  modes  and 
residual  flexibility  modes.  The  residual  flexibility  modes  of  unrestrained  substructures  are 
obtained  from  the  neglected  or  deleted  free  interface  normal  modes  just  like  they  were 
obtained  fit>m  a  substructure  that  is  restrained.  The  only  difference  is  that  the  flexibility 
matrix  obtained  in  Eq.  (26)  is  used.  Residual  flexibility  modes  are  calculated  by  computing 
the  static  flexibility  and  subtracting  the  flexibility  due  to  the  retained  modes.  The  residual 
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flexibility  inodes  are  obtained  fi’om  the  flexibility  matrix  in  Eq.  (26)  in  the  same  way  that 
they  were  obtained  in  Eq.  (25): 


[fj  Ho']-  j  (28) 

It  is  noteworthy  to  state  that  by  performing  the  operation  in  Eq.  (25)  and  Eq.  (28)  that  the 
‘*01”  and  “11”  partitions  of  the  residual  flexibility  matrix  are  extracted  to  form  the  residual 
flexibility  modes. 

C.  BRIEF  HISTORICAL  REVIEW 

The  following  is  a  brief  historical  background  in  the  development  of  componem  mode 
synthesis: 

•  [1965]  In  his  paper.  Dynamic  Amdyas  of  Structund  Systems  by  Componem 
Modes,  Hurty  devdoped  the  first  substructure  coupling  method.  His  technique 
involved  fixed  imeriace  normal  modes,  rigid-body  modes,  and  static  constraint 
modes.  [Ref  5] 

•  [1967]  Bamford  first  introduced  the  concept  of  flexibility  nKxles  in  his  paper  A 
Modal  Combination  Program  for  Dynamic  Analysis  of  Structures.  [Ref  6] 

•  [1968]  In  Cotqtling  of  Substructures  for  Dynamic  Ancdysis,  [Ref  7]  Craig  & 
Hampton  extended  Hurty's  concepts  by  showing  that  rigid-body  modes  did  not 
need  to  be  separated  fium  static  constraim  modes,  but  could  be  calculated  using 
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tlw  same  procedure.  The  Craig-Bampton  procedure  is  oat  of  the  techniques  used 
in  this  study. 

•  [1969]  Goldnuui,  in  his  paper  Vibration  Analysis  of  Dynamic  Partitioning  [Rd!  8] 
and  Hou,  in  his  piq)er  Review  of  Modal  Synthesis  Techniques  and  a  New 
ApprocKh,  [Ref  9]  first  introduced  the  use  of  firee-interface  normal  modes. 

•  [1971]  MacNeal  used  both  free-imerfiice  normal  modes  and  readual  flexibility 
modes  to  couple  substructures.  He  also  suggested  the  use  of  statically  derived 
modes  in  describing  substructure  motion.  These  methods  were  introduced  in  his 
paper  A  Hybrid  Method  of  Component  Mode  Synthesis.  [Ref  10]  The  MacNeal 
method  is  another  technique  that  will  be  analyzed  in  this  study. 

•  [1977]  In  his  Doctoral  Dissertation,  A  General  Procedure  for  Substructure 
Coiqtling  in  Dynamic  Analysis,  [Ref.  11]  Chang,  under  the  guidance  of  Craig, 
showed  how  both  fi^imerfiue  normal  modes  and  residual  flexibility  tiKxles  cmild 
be  employed  to  couple  substructures.  The  Craig/Chang  procedure  is  the  third  and 
final  substructure  coupling  technique  to  be  examined  in  this  study. 

Now  that  the  mode  descriptions  have  been  defined,  and  a  brief  history  of  Axdien  the 
various  methods  of  CMS  were  developed,  the  three  methods  to  be  examined  in  this  study 
can  now  be  derived.  It  is  the  intent  of  Chapter  m  to  describe  the  three  CMS  techniques 
and  how  the  various  mode  sh^}es  are  employed  to  syntherize  substructures  into  a  system. 


m.  COMPONENT  MODE  SYNTHESIS  FORMULATIONS 


The  theory  preseated  herdn  is  takoi  directly  and  exclusively  from  references  [1 1  and  12], 

A.  CRAIG-BAMPTON  FORMULATION 

There  are  three  substructure  coupling  procedures  that  serve  as  potential  candidates  to 
be  used  in  the  mast/antenna  synthesis.  The  Craig-Banq)ton  formulation,  the  Craig-Chang 
residual  flexibility  method,  and  the  MacNeal  residual  flexibility  method.  In  this  section, 
the  Craig-Bampton  formulation  will  be  presorted,  while  in  Section  B  the  Craig-Chang  and 
MacNeal  residual  flexibility  methods  will  be  presented  together  because  of  the  similarities 
in  the  methods. 

The  Craig-Banqnon  reduction  procedure  uses  a  combination  of  static  constraint 
modes  and  fixed  inter&ce  normal  modes  to  reduce  the  conqronent  modd.  Both  the  static 
constraint  modes  and  the  fixed  interfiice  normal  modes  are  obtained  from  the  finite 
dement  substructure  models.  This  combLied  set  of  mode  shapes  will  be  used  to  transform 
the  original  large  order  substructure  mass  and  stiffiiess  matrices  down  to  a  significantly 
smaller  size,  a  size  equd  to  the  number  of  mode  shapes  included  in  the  transformation 
matrix.  The  transformation  matrix  [t,],  for  the  Craig-Bampton  formulation,  contains  the 
shape  fiinctions  as  its  columns  as  follows; 

fe(-P 
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(29) 


This  tnmsfomiation  matrix  is  obtained  for  each  substructure  in  tl»  system.  The  size  of  the 
static  constraint  mode  partition  of  the  transformation  matrix  is  always  fixed  because  the 
number  of  columns  corresponds  to  the  nun:d>er  of  interface  degrms  of  freedom.  Howevo*, 
the  size  of  the  fixed  interfiu:e  normal  mode  partition  is  not  held  constant.  The  number  of 
columns  can  range  as  low  as  one  column  if  only  one  fixed  interface  normal  mode  is 
retained,  or  as  high  "m"  columns  where  "m”  is  the  total  number  of  internal  degrees  of 
freedom.  The  size  of  the  transformation  matrix  depends  upon  how  many  modes  are 
required  to  accurately  represent  the  physical  dynamic  response  of  the  system  when 
subjected  to  a  forced  input.  Retaining  fewer  modes  than  the  total  possible  modes  available 
is  referred  to  as  "modal  truncation,"  and  pro^ades  the  computational  efficioicy  of  the 
method.  Retaining  fewer  modes  than  the  total  amount  of  modes  available  means  fewer 
calculations  required  in  conducting  the  dynamic  analysis.  On  the  other  hand,  if  the  number 
of  modes  retained  are  not  su£5cient  to  accurately  determine  the  dynamic  response,  then 
the  benefits  of  reduced  compute  times  do  not  outwdgh  the  magnitude  of  error  obtained  in 
the  analysis.  Therefore,  while  the  benefits  of  modal  truncation  are  important  in  shortening 
compute  times,  they  are  not  as  important  as  obtaining  accurate  results.  In  terms  of 
computational  efficiency,  large  benefits  can  be  achieved  using  this  method  if  only  the 
lower  range  of  fi^quencies  is  of  int^est.  This  method  is  applicable  to  the  mast  which  is 
subjected  to  typically  low  forcing  frequencies.  By  retaining  a  few  of  each  of  the 
componoit  modes,  an  accurate  assessment  of  the  dynamic  response  of  the  mast  and 
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antenna  is  obtained.  The  examples  contained  in  Chapters  IV  and  V  demonstrate  how  by 
training  ju^  a  few  modes  of  each  substructure,  accurate  results  are  achieved. 

By  pre-muhiplying  the  respective  mass  and  stiffiiess  matrices  by  the  transpose  of  the 
transformation  matiix  and  then  post  multiplying  the  mass  and  stiffiiess  matrices  by  the 
transformation  matrix,  the  reduced  componem  model  is  obtained  as  follows: 


Carrying  out  the  operations  in  Eq.(30a,b)  and  simplifying,  the  partitions  of  the  reduced 
mass  and  stiffiiess  matrices  are  expressed  as  follows: 

k;  -k;'  =(oi 
k;  =[k.  -k„k„-'k„] 
m:  =[i«] 

m;  =m;’  =[»"'(M„(^i,K„,)+M„,)] 

m;  +M„,)  +M„(-Ki,K„,)  +M.] 
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The  tom  'Yeduced”,  designated  by  the  subscript  “r”,  means  that  the  resulting  mass  and 
stiffiiess  matrices,  [K,]  &  [M,],  are  of  smaller  dimension  than  the  original  matrix. 
Although  this  transformation  matrix  reduces  the  size  of  the  component  model,  it  does  not 
assemble  the  individual  substructure  models.  There  is  a  second  transformation  matrix  that 
synthesizes  the  substructures  to  produce  the  total  system  by  enforcing  compatibility  and 
equilibrium  of  intCTfiue  coordinates  and  inter&ce  forces  as  follows; 


WhW) 

(31a) 

(31b) 

The  compatibility  of  interlace  coordinates  denoted  by  Eq.  (31a)  implies  that  the 
displacement  at  the  inter&ce  of  structure  1  equals  the  displacement  at  the  inter&ce  of 
structure  2.  Likewise,  the  equilibrium  of  interface  forces  denoted  by  Eq.  (31b)  implies 
that  the  sum  of  the  forces  at  the  interface  are  equal  but  acting  in  opposite  directions. 

The  generalized  coordinates  employed  in  many  CMS  methods  can  be  identified  with 
the  imerface  DOF  and  the  interior  DOF.  In  order  to  synthesize  two  substructures  into  a 
system,  a  linear  transformation  that  maps  the  set  of  lineariy  independent  coordinates,  p„ 
into  the  set  of  generalized  coordinates  “p”  is  defined  as  follows: 


<p}=fc]4>.} 


(32) 
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Additionally,  the  compatibility  equation,  Eq.  (31a)  can  be  written  in  terms  of  the 
genmlized  coordinates  "‘p”  and  combined  to  form  a  matrix  compatibility  equatio of  the 
form; 


[Cl4)}={0}  (33) 

Let  the  vector  of  generalized  coordinates  be  partitioned  into  linearly  independent  and 
linearly  dependent  coordinates,  and  partition  the  compatibility  matrix  accordingly. 
Equation  (33)  can  now  be  written  as  follows; 

c«]|p‘|=fo)  (34) 

where  subscript  “d”  and  subscript  T’  represent  linear  dependence  and  linear  independence 
respectively.  The  partition  of  the  compatibility  matrix  is  a  nonsingular  square  matrix. 
Expanding  Eq.  (34),  the  dependent  coordinates  can  be  written  in  terms  of  the  independent 
coordinates  as  follows; 


(35) 


Using  Eq.  (35),  the  vector  of  generalized  coordinates  can  now  be  given  by; 


By  satisfying  Eq.  (32),  [T,]  can  be  written  as  follows: 


(37) 

This  tnmsfonnation  matrix  is  employed  to  couple  the  uncoupled  mass  and  stiflfoess 
matrices.  [Ref.  12:p.  472]  The  coupling  procedure  for  the  Craig-Bampton  formulation 
will  now  be  derived.  Let 


P'=g}.  p>=g}  (38.b) 

where  pg  represents  the  generalized  internal  coordinates  and  p,  represents  the  generalized 
interfiu»  coordinates  respectivdy.  Ad^onally,  superscript  ”1"  refers  to  substructure  1 
and  superscript  "2"  refers  to  substructure  2. 

Define  p[  as  the  set  of  dependent  coordinates,  and  let  the  set  of  linearly  independoit 
coordinates  be  expressed  as  follows: 


(39) 


Rewriting  the  compatibility  equation  in  terms  of  the  generalized  coordinates,  Eq.  (31a), 
and  usii%  the  relation  in  Eq.  (38a,b)  the  following  is  obtained: 
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Noting  that  the  set  of  dependent  coordinates  are  designated  by  p[  and  using  Eq.  (33),  Eq. 
(40)  can  be  written  as  follows; 


[[-HIM  [I]  M] 


p] 

Po 

Pl 

Po 


Using  Eq.  (34), 


[c«]=[-«,  E:a.]=[[ol  [I]  [0]] 


(41) 


(42) 


From  Eq.  (42)  and  using  Eq.  (37),  the  second  transformation  matrix  is  obtained  as 
follows: 


fO]  [I]  [0] 
[1]  [0]  [0] 
[0]  [I]  [0] 
[0]  [0]  [I]. 


(43) 


In  order  to  be  consistent  with  the  partitioning  of  the  mass  and  stiffiiess  matrices,  the 
second  transformation  matrix  can  be  rewritten  as  follows; 


pl  [0]  [0] 

ppi  |0]  [I]  [0] 
[0]  [I] 
10]  [I]  [0] 


(44) 
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Equation  (33)  is  now  satisfied  as  follows: 


Po 

in  [0]  [0]' 

pi 

[0]  [I]  [0] 

Po 

1 

Po 

>  s 

[0]  [0]  [I] 

Pi 

_2 

P? 

[0]  [I]  [0]. 

Po  . 

(45) 


The  system  model  is  obtained  fi^om  an  uncoupled  mass  and  stifihess  matrix.  These 
uncoupled  mass  and  stififiiess  matrices  are  thmselves  formed  fi'om  the  reduced  mass  and 
stiffiiess  matrices  fi'om  each  substructure  as  follows: 

m;]  k;] 

where  the  subscript  denotes  “uncoupled.” 

The  coupled  system  mass  and  stiffiiess  matrices  are  obtained  by  pre-multiplying  the 
uncoupled  mass  and  stifiSiess  matrices  in  Eq.  (46a,b)  by  the  transpose  of  the 
transformation  matrix  in  Eq.  (45)  and  then  post-multiplying  the  uncoupled  mass  and 
sdfifiiess  matrices  by  the  transformation  matrix. 


where  the  subscript  “s”  denotes  “system.” 
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Carrying  out  the  op«:ations  in  Eq.  (47a,b)  the  resulting  coupled  system  mass  and  stiffiiess 
matrices  are  given  by; 


[M.]  = 

M’kk 

0 

0 

Mki 

ML 

W= 

K* 

^KK 

0 

0 

KLc 

o  o 

(48a,b) 

M'fK 

Mn. 

0 

0 

K„J 

The  respective  partitions  of  the  system  mass  and  stififhess  matrices  are  expressed  as 
follows: 


^KK  “&!«] 

MJ*  ^Ikk] 

Mic.  +Mi,)] 

ML  +ML)1 

M.  -6-«i,K„,)'^(ML,(-Ki,K<„)’  +M!,,) +ML(-Ki,K„,)'  +Ml,+.. 

(-Ki,K„.)’\Mj„(-Ki,K„)‘  +Mj,)+ML(^i,K<„)'  4Mi] 

Klat 

K.  =^u  -KLiC'Ki.  -KL^i’KL] 

There  are  a  number  of  advantages  to  the  Craig-Bampton  component  mode 
representation.  The  first,  which  is  especially  beneficial  to  the  analysis  of  the  mast/antenna 
system,  is  that  the  reduced  DOF  system  contain  the  interfiice  DOF  explicitly.  This  makes 
it  very  easy  to  couple  mast  and  antenna  substructures.  In  the  figure  on  page  32  is  an 
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illustration  of  the  mast/antenna  system  used  for  the  examples  in  this  thesis.  Along  the 
length  of  tlw  cross  bar  are  various  node  portions.  These  node  positions  serve  to  connect 
the  beams  that  represem  the  cross  bar,  and  can  also  serve  as  nodes  to  connect  antenna  to 
the  crossbar.  By  specifying  different  "connection  coordinates”  (i.e.  the  ”1”  coordinates), 
and  with  separate  mast  and  antenna  data  files  contained  in  the  library,  the  engmea*  can 
quickly  couple  various  antennae  with  tte  mast  and  nq)idly  detennine  the  dynamic 
response. 
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Should  the  location  of  the  antenna  placement  not  be  suitable,  the  o^ineer  can  specify  a 
new  set  of  interface  coordinates  along  the  cross  bar,  plug  in  the  antenna  at  'he  new 
location,  and  rapidly  calculate  a  system  from  which  a  new  dynamic  response  can  be 
calculated. 

Because  of  these  advantages,  and  because  the  Craig-Bampton  component  mode 
representation  tends  to  result  in  accurate  system  frequencies,  as  will  be  shown,  this  is  a 
widely  used  method.  Additionally,  the  NASTRAN  superelement  schone  uses  the  Craig- 
Bampton  component  mode  representation  with  minor  extensions  as  a  solution  path  to  the 
dynamic  response  problem. 

a  CRAIG-CHANG  AND  MACNEAL  RESIDUAL  FLEXIBILITY 

FORMULATIONS 

The  Craig-Chang  and  MacNeal  residual  flexbility  formulations  will  now  be  discussed. 
Due  to  the  similiarity  in  the  methods,  the  Craig-Chang  procedure  will  be  presented  first, 
and  the  modification  of  this  method  to  produce  the  final  system  of  equations  of  the 
MacNeal  method  will  be  discussed  subsequently.  While  the  Craig-Bampton  representation 
uses  a  combination  of  static  constraint  modes  and  fixed  interface  normal  modes,  the  Craig- 
Chang  residual  flexibility  formulation  cond>ines  fi'ee  inter&ce  normal  modes  with  residual 
flexibility  modes;  thus  the  name:  residual  flexibility  method.  The  transformation  matrix 
vdiich  is  used  to  reduce  the  component  mass  and  stiffiiess  matrices  contains  columns  of 
the  retained  or  kept  fiee  interface  normal  modes  and  readual  flexibility  modes. 
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The  tnnsfonnation  matrix,  or  is  shown  as  follows; 


(49) 


In  the  santt  way  that  the  Craig-Bampton  components  were  reduced,  so  to  are  the 
Criig-Chan^  components  reduced  except  now  using  the  transformation  matrix  of  Eq. 
(49).  The  rt...  mass  and  stiffiiess  matrices  are  obtained  as  follows; 


Carrying  out  the  operations  of  Eq.(SOa,b)  the  following  is  obtained  for  the  respective 
partitions  of  the  reduced  mass  and  stifiOiess  matrices; 


m;  =[i„] 

M|  =M|'  =[0„] 

k;  =['a„j 

K;  =Kf  “Bln.] 
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Again,  each  component  has  its  own  transformation  matrix.  A  reduction  in  component 
matrice  size  is  achieved  by  retaining  less  than  the  total  number  of  free  into’fru^e  normal 
modes.  The  number  of  residual  flexibility  mode  ^u^)es  is  fixed,  and  equals  the  number  of 
interface  d^rees  of  fi^om. 

It  is  not  readii;  apparent  how  the  lower  right  hand  partition  of  the  reduced  mass  and 
stiffiiess  matrices  are  obtained.  In  Chapter  II,  the  flexibility  matrix  was  obtained  by 
inverting  the  stiffiiess  matrix  for  a  fully  restrained  substructure.  If  the  substructure  was 
not  restraiiwd,  then  the  flexibility  matrix  was  obtained  by  performing  an  inertia  rdirf 
solution.  However,  it  was  also  shown  that  the  flexibility  matrix  could  be  obtained  by  the 
following  equation: 


(23) 


Additionally,  the  readual  flexibility  nutrix  was  defined  by  the  following  equation: 


The  residual  flexibility  modes  were  obtained  by  post  multiplying  Eq.  (24)  by 


operation  extracts  the  ‘*OF  and  “II”  partitions  of  the  residual  flexibility  matrix. 
Additionally,  the  reridual  flexibility  modes  were  obtained  from  the  free  interface  normal 
modes  in  the  following  way: 


(51) 
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Using  the  relation  between  the  reskhial  flexibility  modes  and  the  deleted  fiiee  interftce 
normal  modes  in  Eq.  (SI),  the  lower  right  partitions  of  the  reduced  mass  and  stiffiiess 
matrices  wiU  now  be  derived. 


Ml  (52) 

Using  the  properties  of  orthogonality: 

(53) 

Substituting  Eq.  (S3)  into  Eq.  (S2)  the  Iowa*  right  partition  of  the  reduced  mass  matrix  is 
obtained  as  follows: 


M;  (54) 

Likewise  for  the  stiflBiess  matrix,  using  Eq.  (S2)  and  the  relationship  between  residual 
flexibility  modes  and  the  deleted  flee  intorftce  normal  modes,  the  following  equation  is 
obtained: 


Because  the  fi'ee  interflice  normal  modes  are  unity  modal  mass  normalized,  the  following  is 
obtained: 


fr:.r[K.H!.]=[A„,.]  (56) 
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Substitutiiig  Eq.  (56)  into  Eq.  (SS),  the  lower  right  partition  of  the  reduced  stiffiiess  matrix 
is  obtained  as  follows: 


k;  (s?) 

Just  like  the  Craig-Bampton  formulation,  each  of  the  substructures  have  their  own 
transformation  matrix.  The  purpose  of  this  transformation  matrix  is  to  reduce  the 
respective  componoit  models  before  synthesis  imo  the  system  modd.  This  transformation 
matrix  does  not  synthesize  the  substructures.  Another  transformation  matrix  is  employed 
to  synthesize  the  substructures.  As  in  tlw  Craig-Bampton  formulation,  this  second 
transformation  matrix  results  from  satisfying  compatibility  and  equilibrium  equations; 


WhW) 

(58a) 

II 

(58b) 

Unlike  the  Craig-Bampton  component  mode  representation  where  the  static  constraint 
modes  are  independent  of  the  fixed  interfrice  normal  modes,  the  residual  flexibility  modes 
are  dependent  upon  the  free  interfrice  normal  modes,  and  a  sinq)le  boolean  matrix  will  not 
synthesize  the  substructures.  This  second  transformation  matrix,  or  [Tj],  will  now  be 
derived. 

The  component  equation  of  motion  in  terms  of  the  phyrical  coordinates  is  given  by 

[M]{x}+{K]{x}={F}  (59) 


37 


Transfonnii^  Eq.  (S9)  to  component  generalized  coordinates  by  letting 


w  “[tKp}  (60) 

the  following  uncoupled  equations  are  obtained; 

[MK.c]feK}-t{Knc]4'K}=t»Sr{F)  (61a) 

(«ll>) 

where  [MKKl[KKKl[^iDo]^<^  the  upper  left  and  lower  right  partitions  of  the 

reduced  mass  and  stifBiess  matrices  respectively,  or  M|  ,K|  ,Mf  and  partitions  of  the 
reduced  mass  and  stififiiess  matrices.  Equations(61a,b)  are  uncoupled  because  residual 
flexibility  modes  are  obtained  as  linear  combinations  of  the  deleted  flee  interface  nonnal 
modes.  The  deleted  flee  interface  normal  modes  are  orthogonal  to  the  kept  free  interface 
normal  modes.  The  response  of  the  ddeted  generalized  coordinates  will  now  be 
approximated  by  the  pseudostatic  response  by  ignoring  the  acceleration  of  the  deleted 
generalized  coordinates  as  follows: 


(62) 


But  it  was  shown  that; 


[K»]  Hk|]  [KnI*.,]  J'Kf  (63) 
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Equations  (62)  and  (63)  can  be  combined  to  give: 


where  the  deleted  generalized  coordinates  are  dependent  upon  the  kept  generalized 
coordinates. 

Combining  Eqs.  (58a),  (66a,b)  and  (67a,b)  with  [T,],  the  compatibility  matrix  in  Eq.  (34) 
is  obtained  as  follows: 


[Cl  =[c„ 


Lw  w  i  (03 


K1] 

[0]  . 


(68) 


From  Eq.  (37), 


(37) 


Since  [C44]  was  defined  as  a  non-singular  square  matrix,  implies  that  [C^]  is  invertible. 
The  inverse  of  [C^^]  is  obtained  as  follows: 


where  [k,]=[i^i>  +i^i)]  '. 
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Combining  Eqs.  (37),  (68)  and  (69)  the  second  transformation  matrix  is  obtained  as 
follows; 


PrU 

[I]  [0] 

[0]  [I]  . 


(70) 


As  in  the  Craig-Bampton  procedure,  the  uncoupled  mass  and  stifBiess  matrices  are 
formed  from  the  reduced  component  mass  and  stifihess  matrices. 


By  pre*  and  post-multiplying  both  the  uncoupled  mass  matrix  represented  by  Eq.  (71a) 
and  the  uncoupled  stifl&iess  matrix  represented  by  Eq.  (71b)  by  [T,]^  and  [Tj] 
respectively,  the  system  equations  of  motion  are  obtained  as  follows; 
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The  partitions  of  the  system  mass  and  stiffiiess  matrices  rq)resemed  by  Eq.  (72a,b)  are 
expressed  as  follows: 


“toe  ^11  “[^\kk 

M,2=M2,  K,2  =K21 

^22  ~0kK  1^22  ’^K^I^Ik] 

where  m,  =k,(M‘DD  +M^)k,.  The  in^a  due  to  high  order  free  interface  normal  modes 
is  represented  by  “m,.”  [Ref.  12:p.  491] 

As  stated  in  the  beginning  of  this  section,  the  MacNeal  component  mode 
representation  would  be  presented.  By  neglecting  the  inertia  due  to  high-order  free 
interface  normal  modes,  (i.e.  “mt”)>  one  obtains  the  MacNeal  mass  and  stiffiiess  matrices 
[Ref  1 1  :p.  59]  as  follows: 

M„  =[I|ck] 

M„  =Ml,  ={0]  K„  =Kl, 

Mm  =[iLc]  K„  =['A,J^  +*5lk,*K] 

As  will  be  demonstrated  through  the  examples,  the  effect  of  neglecting  “m,”  is 
important  when  predicting  the  higher  frequencies.  The  MacNeal  representation  is  accurate 
in  the  lower  and  mid  frequency  range,  but  less  accurate  in  the  higher  frequency  range. 

Just  like  the  Craig-Bampton  component  mode  representation  tho’e  are  several 
advantages  to  both  residual  flexibility  methods.  By  analyzing  both  the  Craig-Chang  and 
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MacNeal  system  mass  and  stiffiiess  matrices,  one  can  see  that  the  final  system  coordinates 
are  just  the  free  interface  normal  mode  coordinates  from  each  substructure.  This  resulted 
from  the  operations  that  were  conducted  in  forming  the  system  mass  and  stiffriess 
matrices.  Additionally,  since  the  residual  flexibility  modes  account  for  the  static  flexibility 
of  all  modes,  the  methods  are  statically  exact.  The  procedure  is  applicable  to  the  mast  and 
antenna  problem.  Since  the  connection  coordinates  are  not  explicitly  retained,  this  allows 
the  engineer  to  retain  fewer  mode  shapes  and  still  obtain  accurate  results.  This  benefit  can 
be  achieved  when  coupling  antenna  systems  such  as  the  SPS*48E  radar  to  the  mast. 
Unlike  a  whip  amenna  which  essentially  is  coupled  at  one  node,  the  SPS-48E  radar  or 
similar  radar  are  connected  to  the  mast  at  more  than  one  node.  Since  there  are  multiple 
connection  points,  and  the  residual  flexibility  method  does  not  retain  the  connection 
coordinates  explicitly,  a  reduction  in  compute  times  result,  an  advantage  not  found  in 
other  substructure  coupling  procedures. 
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IV.  NUMERICAL  VERIFICATION 


In  Chapters  I  through  III,  the  various  types  of  vibrational  mode  shapes  were  defined 
and  used  in  the  derivations  of  the  three  CMS  formulations  that  are  presented  in  this  thesis. 
In  this  chapter,  some  numerical  convergence  examples  are  made  to  compare  system 
natural  frequencies  obtained  by  the  three  CMS  methods  presented  in  this  thesis. 

The  first  example  is  provided  solely  as  a  validation  of  both  the  FE  and  CMS  computer 
codes.  In  his  Doctoral  Dissertation,  Chang  performed  a  numerical  example  by  which  two 
fixed’fm  cantilevered  beams  were  synthesized  into  a  fixed-fixed  cantilevered  beam.  The 
beams  were  synthesized  using  the  various  types  of  CMS  methods.  Example  1  is  a 
r^lication  of  that  example. 

In  example  2,  the  mast  and  antenna  substructures,  illustrated  in  Figure  1,  are 
syntherized  imo  the  mast  and  antenna  syston.  The  system  is  synthesized  three  difierent 
times  using  a  varying  amount  of  substructure  mode  shq)es.  The  resulting  natural 
fioquencies  are  obtained  using  MATLAB'S  dgensolver.  A  conq>arison  is  made  between 
each  of  the  methods  and  against  a  FE  generated  modd  of  the  mast  and  antenna  system. 
This  comparison  reflects  both  the  accuracy  of  the  three  methods  in  predicting  the  natural 
fi^uencies  and  the  number  of  calculations,  or  floating  point  operations  "FLOPS",  that  are 
required  to  synthesize  the  structures  and  poform  the  eigensolution. 
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A.  NATURAL  FREQUENCY  CALCULATION  AND  COMPARISON  OF  A 
FIXED-FIXED  CANTILEVERED  BEAM  FOR  PURPOSES  OF  FINITE 
ELEMENT  AND  COMPONENT  MODE  SYNTHESIS  CODE  VALIDATTON. 

As  stated  in  the  beginning  of  this  chapter,  the  purpose  of  this  example  is  to  validate  the 

FE  and  CMS  code  which  were  written  by  the  author.  Figure  2  illustrates  a  fixed-fixed 

cantilevered  beam  that  was  synthesized  fix)m  two  fixed-fiee  cantilevered  beams. 


y\ 

y  ■  --  ^ 

substructure  1 


i  N 

!  ■  ■  X 

I 

I 

substructure  2  \ 


Figure!:  Two  Cantilevered  Beam  Substructures 

Each  substructure  has  two  degrees  of  freedom  at  each  node;  one  translational  and  one 
rotational.  The  first  substructure  is  represented  by  4  beam  elements  corresponding  to  8 
DOF,  while  the  second  substructure  is  represrated  by  3  beam  elemems  corresponding  to  6 
DOF.  After  synthesis,  the  system  has  a  total  of  12  DOF.  As  in  Chang's  example,  each 
method  was  run  using  a  total  of  4,  6,  8  and  10  DOF,  or  33,  50,  67,  and  83%  of  all 
available  modes.  The  results  contained  in  Tables  1-4  are  expressed  in  "Hz"  and  are  the 
exact  same  results  that  are  found  in  Tables  12-15  of  the  referenced  example  ^^ch  are 
expressed  in  “(rad/sec)^ .” 

In  analyzing  the  results  contained  in  Tables  1-4,  it  spears,  for  this  case,  that  the 
Craig-Chang  procedure  yielded  a  better  prediction  of  the  natural  frequencies  than  the 
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other  mo  fbmnilations.  The  results  are  summarized  in  Figure  3:  Comparison  of  CMS 
methods  with  fi^uency  (Hz)  error  of  <=  0. 1  percent  for  a  fixed-fixed  beam.  In  all  cases, 
the  MacNeal  procedure  predicted  a  natural  frequency  in  access  of  fifty  percent  error  as 
compared  to  the  standard  finite  element  solution.  This  error  in  the  predicted  natural 
frequency  occurred  when  predicting  the  highest  mode  only.  In  general,  the  MacNeal 
procedure  predicted  the  lower  range  of  natural  frequencies  comparativdy  well,  but 
appeared  to  faher  when  predicting  the  higl^  frequency  range.  A  possible  reason  for  this 
occurence  comes  as  a  result  of  n^ecting  the  inertia  of  high  order  free  inter&ce  normal 
modes,  or  "mj"  as  delineated  in  Chapter  HI. 
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TABLE  3.  COUPLED  SYSTEM  NATURAL  FREQUENCY  COMPARISON  (FIXED-FIXED  BEAM) 

Percent  of  Substnictiire  1  Modes  R^nined:  63% 
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Fignre  3:  Comparison  of  CMS  Methods  With  Frequency  (Hz) 
Error  of  <=  0.1  Percent  for  ‘^Fixed-Fixed’*  Beam  System 
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B.  NATURAL  FREQUENCY  CALCULATION  AND  COMPARISON  OF 

MAST  AND  ANTENNA  SYSTEM  MODEL 

The  purpose  of  the  first  example  was  to  validate  both  the  FE  code  and  the  various 
CMS  codes.  The  CMS  code,  written  in  MATLAB,  will  accommodate  the  most  general 
FE  modeling.  In  other  words,  the  CMS  code  is  written  indep^oit  of  the  type  of  FE 
modeling.  Should  the  mast  and  antenna  system  be  modeled  by  beam  elements  with  6 
DOF/node,  tlw  CMS  code  can  be  q)plied  to  it  without  any  modification. 

In  this  section  another  numerical  example  which  demonstrates  the  use  of  tte  three 
component  mode  synthesis  procedures  is  presemed.  A  "mock-up"  mast  and  antenna 
system  consisting  of  17  dements  was  assonbled  using  a  standard  finite  demem  procedure 
and  aU  three  component  mode  synthesis  procedures.  Although  the  mast/antenna  system 
that  is  used  in  the  exanq;>les  is  modded  with  just  a  few  dements,  the  resulting  models  are 
laii^  enough  to  allow  the  effects  of  mode  truncation  to  be  assessed.  Again,  it  is  not  the 
imem  of  this  report  to  solve  a  base  exdtadon  problem  on  a  realistic  mast  and  amenna 
model,  but  rather  to  dmnonstrate  how  CMS  can  be  used  when  performing  dynamic 
analyses  for  design  purposes. 

Just  like  the  first  exan^le,  this  example  compares  natural  fi^uency  calculations  for 
the  total  mast/amenna  system  as  computed  using  the  three  CMS  formulations  as  wdl  as 
using  a  standard  FE  procedure.  Tables  5-7  contain  the  results  of  this  comparison.  In  the 
tables,  each  row  contains  the  estimate  of  a  mode  fi^uency.  The  first  column  contains  the 
mast/amenna  system  natural  frequency  estimates  as  calculated  using  the  standard  FE 
procedure,  and  serves  as  the  reference  value  against  which  the  CMS  natural  fiequency 


52 


estimates  are  to  be  compared.  Columns  2  through  4  contain  the  analogous  natural 
fiequency  estimates  and  percent  error,  as  calculated  from  the  nuist/antenna  systei  i  model 
synthesized  using  each  of  the  three  CMS  procedures.  Also  included  in  the  column 
headings  are  floating  point  operations  (FLOPS)  coums  which  provide  a  measure  of  the 
number  of  calculations  required  to  assemble  the  mast/antenna  syston  and  calculate  the 
natural  frequencies  and  mode  sluq)es. 

Table  1  presoits  the  system  frequency  comparison  where  18%  of  the  available  mast 
modes  are  retained,  and  22%  of  the  available  amoina  modes  are  retained.  Table  2  repeats 
the  calculations  with  42%  of  the  available  mast  modes  retained  and  39%  of  the  available 
antenna  modes  retained,  and  Table  3  repeats  the  calculations  with  79%  of  the  available 
mast  modes  retained  and  67%  of  the  available  antenna  modes  retained.  Note  that  each 
subsequent  table  presents  comparisons  for  an  increasing  number  of  mode  frequencies  due 
to  the  fact  that  an  increase  in  the  number  of  retained  componem  modes  makes  possible  an 
increase  in  the  number  of  system  modes  which  may  be  calculated. 

Note  that  in  the  FEM  model,  the  FLOPS  coum  stays  fixed  at  slightly  over  2‘10‘. 
This  is  a  rather  small  number  as  the  model  is  a  small  model  when  compared  to  one  that  a 
design  engineer  would  generate  for  analysis  of  an  actual  mast/antenna  assembly.  It  is 
notewortl^  to  state  that  in  tlus  particular  model,  0.7 ’10‘  FLOPS  were  expended  in 
computing  the  combination  of  fixed  interface  normal  modes  and  static  constraint  modes 
u^g  the  Craig-Bampton  procedure.  Additionally,  0.8  *10*  FLOPS  were  expended  in 
computing  the  free  interface  normal  modes  uting  the  Craig-Chang  and  MacNeal 
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procedures.  Theoretically,  once  the  various  vibrational  modes  have  been  found,  they  med 
not  be  calculated  again.  Note  also  that  in  all  three  models,  these  figures  comprise  a 
significant  portion  of  the  total  FLOPS. 

From  Tables  S  through  7,  it  is  seen  that  all  three  methods  produce  excellent  frequency 
predictions.  All  three  methods  demonstrate  sudden  increases  in  frequency  error  above  a 
certain  mode.  This  reveals  the  extent  to  which  the  retained  componrat  modes  accurately 
represent  the  dynamics  of  the  synthesized  mast/antenna  system.  In  the  MacNeal 
procedure,  the  percentage  error  exceeded  100%  when  calculating  the  highest  mode.  This 
error  in  predicting  the  highest  frequency  mode  could  possibly  be  attributed  to  neglecting 
the  inertia  due  to  high-order  free  interface  normal  modes.  Note  that  by  neglecting  the 
inertia,  that  the  accuracy  in  predicting  natural  fiYX}uencies  is  only  effected  at  the  last  few 
modes. 

As  stated  in  the  previous  paragraph,  ail  three  methods  produced  excellent  results  in 
predicting  natural  fiequencies,  but  with  less  cost  in  terms  of  number  of  computations  as 
compared  to  the  standard  FE  calculation.  The  Craig-Chang  procedure  in  general  provided 
the  greatest  number  of  natural  fi^uencies  with  error  less  than  or  equal  to  0.1%  (in  Hz) 
(see  Figure  4).  However,  the  Craig-Bampton  procedure  yielded  the  same  number  of 
fiequencies  with  error  less  than  or  equal  to  0.1%  as  the  Craig-Chang  procedure  whoi 
retaining  a  laige  of  number  of  component  modes,  but  at  a  slightly  more  cost  than  the 
Craig-Chang  procedure. 
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TABLE  5.  COUPLED  SYSTEM  NATURAL  FREQUENCY  COMPARISON  (MAST/ANTENNA  MODEL) 

Percent  of  Available  Mast  Modes  Retained:  18% 

Percent  of  Available  Antenna  Modes  Retained:  22% _ 
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TABLE  6.  COUPLED  SYSTEM  NATURAL  FREQUENCY  COMPARISON  (MAST/ANTENNA  MODEL) 

Percent  of  Available  Mast  Modes  Retained:  42% 

Percent  of  Available  Antenna  Modes  Retained:  39% _ 
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Notes:  The  column  “Standard  Finite  Element”  contains  the  mode  frequencies  for  the  mast/antenna  system  assembled  using  a 
standard  finite  element  approach.  The  corresponding  mode  frequency  estimates  for  the  three  synthesis  methods  are  shown  in 
the  three  colunms  to  the  right.  The  Floating  Point  Operations  (n..0PS)  required  for  each  synthesis  method  and  the  finite 
element  calculation  are  shown.  Retained  mode  percentages  are  rounded  to  the  nearest  integer  for  clarity.  Modes  1 1  through 
18  excluded  to  conserve  space. 


TABLE  7.  COUPLED  SYSTEM  NATURAL  FREQUENCY  COMPARISON  (MAST/ANTENNA  MODEL) 

Percent  of  Available  Mast  Modes  Retained:  79% 

Percent  of  Available  Antenna  Modes  Retained:  67% _ _ 
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Notes;  The  column  “Standard  Finite  Element”  contains  the  mode  frequencies  for  the  mast/antenna  system  assembled  using  a 
standard  finite  element  approach.  The  corresponding  mode  frequency  estimates  for  the  three  synthesis  methods  are  shown  in 
the  three  columns  to  the  right.  The  Floating  Point  Operations  (n..OPS)  required  for  each  synthesis  method  and  the  finite 
element  calculation  are  shown.  Retained  mode  percentages  are  rounded  to  the  nearest  integer  for  clarity.  Modes  6  through  27 
excluded  to  conserve  space. 


0%  10%  20%  30%  40%  S0%  00%  70%  80% 

Reront  of  avaifadile  ^ystea  nodes  Rtaiacd 

Figure  4:  Comparison  of  CMS  Methods  With  Frequency  (Hz) 
Error  of  0.1  Percent  for  Mast/Antenna  System 
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V.  BASE  EXCITATION  FORMULATIONS 


In  the  previous  section,  the  natural  frequencies  of  the  mast  and  antenna  system  were 
calculated  for  increasing  number  of  retained  component  modes.  Natural  frequencies  and 
mode  shapes  are  important  modal  parameters  and  are  fundamental  in  solving  for  the 
forced  re^nse  of  a  system.  As  demonstrated  in  the  previous  section,  accurate  natural 
frequencies  of  a  system  can  be  obtained  using  CMS  at  a  cost  less  than  that  associated  with 
standard  FE  modeling. 

In  this  section,  two  base  excitation  formulations  will  be  presented.  The  first 
formulatitMi  requires  the  knowledge  of  the  acceleration  of  the  mast  base  coordinates  (i.e. 
the  coordinates  where  the  mast  and  ship  are  coupled)  as  a  function  of  time.  In  other 
words,  the  formulation  requires  that  the  acceleration  time  history  of  the  base  coordinates 
be  known.  The  second  formulation  requires  the  knowledge  of  the  displacement  of  the 
base  coordinates  as  a  function  of  time,  or  the  displacemem  time  history  of  the  base 
coordinates.  Using  both  formulations,  nummical  convergence  assessments  will  be  made, 
and  the  benefits  that  CMS  has  to  ofiEer  the  rnast/amenna  design  process  will  be 
demonstrated. 
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A.  BASE  EXCITATION  FROM  PRESCRIBED  ACCELERATION 

Once  the  FE  program  has  numerically  assembled  the  mast  mass  and  stiffiiess  matrices, 
and  the  acceleration  of  the  base  coordinates  are  specified  as  a  function  of  time,  the  base 
excitation  pr(ri)lem  can  be  derived  fiom  the  following  equation  of  motion: 


where  subscript  "O"  represents  the  interior  coordinates,  and  subscript  "B"  represents  the 
base  coordinates. 

Solving  the  top  row  of  equations  in  Eq.  (73)  the  following  is  obtained: 

^00*0 

Since  the  acceleration  of  the  base  is  presoibed,  the  base  acceleration  term  will  be  moved 
to  the  right  hand  of  the  equals  sign  to  obtain  the  following: 

h^OO^O  ~^0  ^OB^B  (^5) 

From  the  bottom  row,  the  fi^Uowing  equation  is  obtained: 

^BO^O  “^^BB^B  "^I^BO^B  “^I^BB^B 

From  Eq.  (76),  the  following  relation  is  obtained  for  the  base  displacment: 
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Equatkm  (77)  is  now  substituted  into  Eq.  (75),  and  after  simplifyii^  the  following 
eqpiaticMi  of  motion  in  terms  of  the  interior  coordinates  is  given  by; 


+[Ko„  =F„  +JC„Ki,M„  (78) 

Equation  (78)  is  the  system  equation  of  motion  of  the  internal  coordinates  in  terms  of  the 
prescribed  base  acceleration  for  the  mast  only.  Since  the  antenna  has  no  prescribed  base 
acceleration,  the  above  formulation  is  not  applied  to  the  FE  generated  mass  and  stiffiiess 
matrices  of  the  antenna. 

Now  that  both  sid>structures  have  been  numerically  constructed,  they  are  now  ready  to 
be  synthesized  into  the  mast  and  antenna  system.  Both  the  mast  and  antenna  substructure 
mass  and  stHBicM  matrices  and  force  vector  are  now  partitioned  into  internal  and 
connection  coordinates  as  follows: 


\^iere  the  subscript  “O”  and  subscript  “I”  represent  the  internal  and  interface  coordinates 
respectively. 

Using  the  procedures  of  Ch^er  lH,  the  respective  substructure  mass  and  stiffiiess 
matrices  and  force  vectors  are  transformed  into  a  “reduced”  system  by  using  the  first 
transformation  matrix  as  follows; 

[Mj-frJCMJfr,]  [K,]=[r,r(Kl[T,]  ^,}-[r,r{F}  (79a,b,c) 
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As  derived  in  Chapter  III,  the  coupled  mast/antenna  system  mass  and  stiffiiess  matrices 
and  force  vectors  are  formed  from  an  uncoupled  system  uring  the  second  transformation 
matrix,  and  is  given  by: 

Kl=[Tj[M.Ir.]  [K.]-[rJpC.Ir,]  (80a,b,c) 

where  the  subscripts  “s”  and  **u”  represoit  the  coupled  system  and  uncoupled  system 
respectivdy. 

The  operations  in  Eq.  (79  a,  b,  c)  and  Eq.  (80  a,  b,  c)  are  applicable  to  each  of  the 
CMS  methods  using  each  method’s  respective  transformation  matrices,  (i.e.  [r,]  and 
Now  that  the  system  force  vector  and  system  mass  and  sdffiiess  matrices  have  been 
formed,  the  coupled  ^em  natural  frequoides  and  mode  shapes  can  be  obtained  from  the 
following: 

(81) 

In  Chapter  III  it  was  shown  how  the  generalized  coordinates  or,  “p”  coordinates,  were 
obtained  from  a  linear  transformation  of  the  phyrical  coordinates  or,  ‘V’  coordinates, 
using  the  foat  transformation  matrix.  Additionally,  a  set  of  linearly  independent 
coordinates  were  obtained  from  the  set  of  generalized  coordinates,  which  consist  of 

I 

linearly  independent  and  linearly  depends  coordinates,  using  the  second  transformation 
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matrix.  Therefore,  the  coordinates  associated  with  Eq.  (80  a,  b,  c)  and  Eq.  (81)  are  the 
set  of  generalized  linearly  independent  coordinates  or,  “p,  ,”  as  defined  in  Chapt^  III. 

The  set  of  generalize  linearly  independent  coordinates  are  transformed  into  a  set  of  modal 
coordinates  by  the  following  relation: 


{(>o}=[*"]llo)  (K) 

where  }  represents  the  vector  of  generalized  linearly  independent  internal  coordinates, 
^o}  represents  the  vector  of  modal  internal  coordinates,  and  represents  the  matrix 
of  unity  modal  mass  fiiee  interfiux  normal  modes  obtained  fi^m  the  solution  to  Eq.(81). 

The  system  equation  of  motion  is  then  pre-multiplied  by  the  transpose  of  the  matrix  of 
normal  modes  to  obtain  the  following  modal  system  of  equations: 

Wfeo  l+fA.Jfio  }=  {F}  (83) 

}  represents  the  vector  of  modal  accelerations  and  {F}  represents  the  vector  of 
modal  forces. 

There  are  “m”  equations  associated  with  the  solution  to  Eq.  (83)  where  “m”  is  equal 
to  the  number  of  rows  or  columns  in  either  the  system  mass  or  stifihess  matrice.  The 
acceleration  at  the  base  is  a  prescribed  harmonic  input. 
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The  base  acceleration  is  taken  as: 


^,}  =  -0‘^,}sin(0t)  (84) 

where  “Q ”  represents  the  forcing  fiequency  expressed  in  rad/sec,  and  {Xg}  re{Hesents  the 
vector  of  amplitudes  of  the  base  displacement  expressed  in  inches.  Since  accderation  is 
prescribed  at  the  base,  and  the  product  of  mass  and  acceleration  is  consistoit  with  the 
units  of  force,  the  resulting  forcing  function  is  also  harmonic. 

Therefore,  the  form  of  the  ‘1^”  equation  of  this  system  of  equations,  (i.e.  Eq.  (83)),  is 
given  by: 


qo,  =F,(sin(Qt))  (85) 

where  =  X  ,  and  is  associated  with  the  **i^”  natural  frequency  of  the  system. 

Equation  (85)  can  be  solved  u^g  two  different  methods.  Since  Eq.  (85)  is  an 
ordinary  differential  equation,  the  solution  can  be  obtained  by  finding  the  particular 
solution  which  requires  the  knowledge  of  the  initial  conditions,  (i.e.  <l>^  ■q,^(t  =  0)  and 
^  'qio(t  ~  0))  or  the  initial  displacement  and  vdocity  of  the  system.  Since  there  are  no 
initial  conditions  on  the  mast  and  antenna  prior  to  the  shock  wave  meaning  that  the  mast 
and  amenna  is  motionless  prior  to  the  shock  wave,  a  more  convenient  and  preferred  way 
to  solve  Eq.  (85)  is  through  convolution.  The  homogenous  solution  to  Eq.  (85)  is  given 
by: 

qi(t)  =  A,(cos(«it))  +  Bi(sin(«,t))  (86) 
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wiwre  “A,  ”  and  ”  are  the  amplitudes  of  the  iiKxlal  displacement  and  are  evaluated 
from  the  initial  conditions. 

Using  convolution,  the  particular  solution  is  given 

qi(t)  *  A,(cos(a),t))  +  Bj(sin(a},t))  +  f  —  Mn(Wi(t  -t  ))Fjan(0  r  )dT  (87) 
where  “r  “  is  a  time  constant. 

The  im^ral  in  Eq.  (87)  is  generally  knovm  as  the  convolution  integral  or  more  formally  as 
Duhamel’s  int^ral.  After  solving  Duhamd*s  integral,  Eq.  (87)  can  be  rewritten  as 
follows: 


FXt  sin((i>,t)  /^sin(Qt) 

,.(,).A,(cos(«..))^B,(sin(a,..))  - 


(88) 


Since  there  is  no  initial  di^lacement  or  velodty,  “A,  ”  and  “B,  ”  are  both  equal  to  0. 


Equation  (88)  is  then  simplified  as  ftrllows: 


q,(t)=  - 


F,Q  sin(&>tt)  /^sin(Ot) 

(■}j((i)|^  ”0*)  (wi*  “^*) 


(89) 


Equation  (89)  is  solved  for  each  of  the  *‘m”  equations  in  the  system  of  equations.  By 
solving  Eq.  (89),  the  modal  displacemem  for  each  coordinate  is  obtained  as  a  function  of 
time. 

Since  the  mast  was  modeled  as  a  ‘Tree-free”  structure  meaning  that  there  were  no 
restrained  coordinates,  the  synthesis  of  the  mast  and  antenna  system  yields  a  structure  that 
is  also  “fiiee-firee”.  This  means  that  there  are  some  rig^d  body  modes  associated  with  the 
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dgensohition  of  the  system  of  equations.  Since  the  mast  was  modeled  with  three 
DOF/node,  there  are  three  rigid  body  modes  v^diich  correspond  to  three  natural 
frequencies  of  0  rad/sec.  Therefore,  for  the  first  three  equations  of  Eq.  (83),  Eq.  (89) 
cannot  be  applied  in  solvii^  for  the  modal  (tisplacement  of  Xht  first  three  coordinates  as  a 
function  of  tin».  Using  Eq.  (85)  and  substituting  “0”  for  u)  Eq.  (85)  is  simplified  as 
follows: 


qo.(t)=i^sin(0t).  i=l:3  (90) 

Integrating  Eq.  (90),  the  following  is  obtained  for  the  modal  velocity  as  a  fimction  of  time 
for  the  first  three  modal  coordinates: 

')o,(t)”^oo»(Ot)+C|,  i-l:3  (91) 

where  “c”  is  a  constant  of  imegration. 

Integrating  Eq.  (91),  the  following  is  obtained  for  the  modal  displacem^it  as  a  function  of 
time  for  the  first  three  modal  coordinates: 


^o,  (0  =^sin(Ot)  +  Cjt  +  dj,  i=l  :3  (92) 

where  “d”  is  a  second  constant  of  integration. 

The  constants  of  integration  “c”  and  “d”  are  obtained  by  substituting  the  initial  conditions 
into  Eq.  (91)  and  Eq.  (92)  and  solving  Eq.  (91)  and  Eq.  (92)  for  “c”  and  “d”.  Substituting 
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the  initial  condition  of  displacement  into  Eq.  (92)  the  second  constant  of  integration  is 
obtained: 

d,  =0.  i=l:3 

The  first  constant  of  int^ration,  *V*  is  obtained  in  a  similar  manner.  Substitutii^  the 
vdodty  initial  conditions  in  Eq.  (91),  the  first  time  constant  is  <4)tained: 


Substituting  the  constants  of  int^ration  into  Eq.(92),  the  modal  displacements  for  the  first 
three  cotudinates  can  be  obtained  fiom  the  following  relation: 

q^(t)  =^silKOt)+^t.  i-l:3  (93) 

Equation  (93)  is  used  to  obtain  the  modal  displacement  solution  as  a  function  of  time  for 
the  first  three  coordinates,  while  Eq.  (89)  is  used  to  obtain  the  modal  displacement 
solution  as  a  function  of  time  for  the  remainder  of  the  set  of  coordinates. 

Once  the  modal  re^nse  is  obtained,  Eq.  (82)  is  used  to  obtain  the  response  in  linearfy 
indepoident  generalized  coordinates.  By  using  the  second  transformation  matrix  followed 
in  succesaon  by  the  first  transformation  matrix  in  a  manner  similar  to  that  of  Eq.  (82),  tl» 
re^nse  of  all  of  the  physical  coordinates  that  defines  the  system  is  obtained. 
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In  what  follows,  the  relation  in  Eq.  (78)  is  used  in  conjunction  with  a  standard  FE 
model  of  the  total  mast/antenna  system  to  perform  the  prescribed  base  acceleration 
dynamic  analysis.  This  analysis  serves  as  the  reference  against  which  the  results  of  various 
CMS  formulations  are  to  be  compared.  To  restate,  in  each  CMS  formulation,  Eq.  (78)  is 
used  to  define  the  mast  component  model,  as  only  the  mast  has  prescribed  base 
accelerations. 

In  the  following  section,  two  numerical  examples  are  provided.  It  is  the  intent  of  the 
examples  to 

( 1 )  compare  the  results  obtained  from  the  three  CMS  formulations 

(2)  and  demonstrate  the  boiefits  of  u^g  CMS  versus  standard  FE  modeling 
when  solving  base  excitation  problems. 

1.  Tip  Deflection  Calculation 

On  page  70  is  a  diagram  of  the  mast  and  antenna  system  being  subjected  to  base 
excitation.  The  excitation  was  performed  at  two  different  and  arbitrarily  selected 
frequencies.  The  first  frequency  was  at  8.95  Hz.  This  frequency  falls  between  mode  1  and 
mode  2  of  the  total  mast/antenna  system.  The  second  frequency  was  between  mode  9  and 
mode  10  at  219  Hz.  The  wide  spread  in  the  frequencies  was  intended  to  demonstrate  that 
many  more  modes  need  to  be  retained  when  calculating  the  response  to  higher  frequency 
excitation  as  compared  with  lower  frequency  excitations.  In  this  example,  the  antenna  tip 
d^ection  was  calculated  using  the  standard  FE  procedure  and  the  three  CMS  procedures. 
The  percent  error  in  antenna  tip  deflection  was  plotted  versus  the  percent  of  available 
ccmiponent  modes  retained  (see  Figures  6-9).  The  calculations  were  performed  twice.  In 


68 


the  first  calculation,  mast  modes  were  tnmcated  while  retaining  all  of  the  available  antenna 
modes,  and  antenna  modes  were  tnmcated  in  the  second  calculation  while  retain' ng  all  of 
the  available  mast  modes.  When  the  mast  was  subjected  to  the  forced  input  at  the  lower 
fiequency  the  Craig-Bampton  and  Craig-Chang  procedure  yielded  results  which 
converged  more  nq)idly  to  the  exact  answer  than  the  MacNeal  procedure. 
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Although,  hard  to  determine  from  Figures  6  and  7,  the  Craig-Bampton  procedure  yielded 
the  best  results  using  fewer  modes  than  the  Craig-Chang  procedure  in  both  the  mast 
truncation  and  antenna  truncation  runs.  When  20%  of  the  available  mast  modes  were 
used,  all  three  procedures  predicted  a  tip  deflection  measurement  that  was  within  0.0S% 
of  the  "exact”  value  (the  “exact”  value  was  calculated  using  the  standard  FE  procedure). 
This  accurate  assessment  was  obtained  at  a  cost  of  60%  of  the  number  of  calculations  that 
were  required  of  the  FE  solution.  The  error  obtained  when  20%  of  the  available  mast 
modes  were  retained  was  significantly  small.  Since  all  of  the  vibrational  mode  shapes  were 
obtained  in  this  calculation  (i.e.  fixed  interface  normal  modes,  free  interface  normal  modes, 
static  constraint  modes,  and  residual  flexibility  modes),  future  assessments  which  retain 
more  mast  modes  would  come  at  an  even  lesser  cost  than  the  initial  assessment. 

When  the  mast  was  subjected  to  the  higher  forcing  frequency,  all  three  methods 
converged  more  slowly  as  compared  to  the  lower  forcing  frequency  when  tnmcating  both 
mast  modes  and  antenna  modes  (see  Figures  8  and  9).  Since  the  forcing  frequency  was 
higher,  more  modes  needed  to  be  retained  in  order  to  obtain  accurate  results.  From  the 
results  of  the  tip  deflection  calculations,  it  appears  that  the  combination  of  fixed  inter&ce 
normal  modes  and  static  constraint  modes  have  led  to  the  higher  rate  of  convergence  using 
the  Craig-Bampton  procedure.  However,  the  results  obtained  using  the  Craig-Chang 
procedure  compared  quite  well  with  the  results  obtained  using  the  Craig  Bampton  method. 
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Figure  6:  Percent  Error  in  Antrana  Tip  Deflection  Plotted  Versus  the 
Percent  of  Available  Mast  Modes  Retained.  (Forcing  Frequency:  8.95  Hz) 
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Figure  7:  Percent  Error  in  Antenna  Tip  Deflection  Plotted  Versus  the 
Percent  of  Available  Antenna  Modes  Retained.  (Forcing  Frequency:  8.95  Hz) 
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Figure  8:  Percent  Error  in  Antenna  Tip  Deflection  Plotted  Versus  the 
Percent  of  Available  Mast  Modes  Retained.  (Forcing  Frequency:  219  Hz) 
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Figure  9:  Percent  Error  in  Antenna  Tip  Deflection  Plotted  Versus  the 
Percent  of  Available  Antenna  Modes  Retained.  (Forcing  Frequency:  219  Hz) 
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2.  Moment  and  Shear  CalcuhiHon 

In  order  to  assess  structural  survivability,  accurate  prediction  of  the  internal 
stresses  in  the  antenna  and  the  antenna/nuist  inter&ce  must  be  calculated.  Therefore,  this 
section  demonstrates  the  calculation  of  the  internal  peak  dynamic  bending  moments  and 
shear  loads  in  the  antenna. 

The  moment  and  shear  calculations  were  calculated  using  the  three  CMS 
procedures  and  the  results  are  compared  in  the  figures.  The  percent  error  in  moment  and 
shear  were  plotted  versus  the  percent  of  a>aulable  mast  modes  retained  and  perc^  of 
available  antenna  modes  retained.  Again,  the  same  two  forcing  firequencies  used  in  the  tip 
deflection  calculation  are  used  here  in  the  moment  and  shear  calculations;  specifically  8.95 
and  219  Hz. 

The  results  obtained  when  calculating  the  shear  and  moment  at  the  mast/antenna 
connection  mirror  the  results  of  the  tip  deflection  calculations  (see  Figures  10-17).  Again, 
the  Craig-Bampton  procedure  yielded  resuhs  that  converged  more  quickly  to  the  "exact" 
answer  (provided  by  standard  FE  calculations)  than  the  other  methods.  However,  the 
results  obtained  using  the  Craig-Chang  procedure  were  quite  similar  to  those  obtained 
using  the  Craig-Bampton  method.  Despite  a  large  initial  error  produced  by  the  MacNeal 
method  as  compared  to  the  other  two  methods,  neady  "exact”  solutions  were  obtained  at 
a  cost  much  less  than  using  standard  FE  calculation  procedures.  If  the  moment  and  shear 
at  the  mast/antenna  interftce  exceeded  an  ^propriate  failure  criteria,  the  antenna  can  be 
easily  rdocated  fi:om  the  end  node  to  anotho-  node  along  the  cross  bar  by  redefining  the 
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connection  coordinates  of  the  mast.  New  moment  and  shear  calculations  would  be  made 
until  an  acceptable  response  obtained.  Redefining  the  connection  coordinates,  synthesizing 
the  new  structure,  and  calculating  the  response  is  much  more  convenient  and 
computationally  efficient  than  reassembling  the  mast  and  antenna  system,  which  would  be 
required  using  standard  FE  procedures. 

Whoi  the  mast  and  antenna  system  were  subjected  to  the  forced  input  at  the  higher 
forcing  frequency,  the  rate  of  convergence  was  again  much  slower  than  that  which  was 
obtained  at  the  lower  excitation  frequency.  However,  all  methods  yield  accurate  results 
at  a  computational  cost  less  than  using  the  standard  FE  procedure  with  the  higher  forcing 
frequency. 
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Figure  10:  Percent  Error  is  Shear  at  Mast  and  Antenna  Connection  Plotted 
Versus  the  Percent  of  Availahle  Mast  Modes  Retained. 

(Forcing  Frequency:  8.95Hz) 
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Figure  11:  Percent  Error  in  Moment  at  Mast  and  Antenna  Connection  Plotted 

Versus  the  Percent  of  Available  Mast  Modes  Retained. 

(Forcing  Frequency:  8.95  Hz) 
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Figure  12: 


Percent  Error  in  Shear  at  Mast  and  Antenna  Connection  Plotted 
Versus  the  Percent  of  Available  Antenna  Modes  Retained. 
(Forcing  Frequency:  8.9S  Hz) 
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Figure  13:  Percent  Error  in  Moment  at  Mast  and  Antenna  Connection  Plotted 

Versus  the  Percent  of  Available  Antenna  Modes  Retained. 

(Forcing  Frequency:  8.95  Hz) 
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F^re  15:  Percent  Error  in  Moment  at  Mast  and  Antenna  Connection  Plotted 

Versus  The  Percent  of  AvaflaUe  Mast  Modes  Retained. 

(Forcing  Frequency:  219  Hz) 


78 


Percent  error  in  nomeBt  at  mast  Percent  error  in  shear  at  mast  and 

and  antenna  connection  antenna  connectioB 


Percent  error  in  shear  at  asast  and  antenna  connectioB  plotted  versus  the  percent  of 
avaiahle  antenna  arodes  retained 


Number  of  antenna  modes  retained  (in  percent) 


Figure  16:  Percent  Error  in  Shear  at  Mast  and  Antenna  Connection  Plotted 
Versus  the  Percent  of  Available  Antenna  Modes  Retained. 
(Forcing  Frequency:  219  Hz) 
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F^re  17:  Pwcent  Error  in  Moment  at  Mast  and  Antenna  Connection  Plotted 

Versus  the  Percent  of  Available  Antenna  Modes  Retained. 

(Forcing  Frequency:  219  Hz) 
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a  BASE  EXCITATION  FROM  PRESCRIBED  DISPLACEMENT 

The  following  derivation  is  applicable  to  the  base  excitation  problem,  where  bas : 
displacements  (as  opposed  to  base  accelerations)  are  to  be  prescribed. .  As  in  the  base 
excitation  from  prescribed  acceleration,  the  derivation  starts  with  the  FEM  generated  mass 
and  stiffiiess  matrices  as  follows; 


The  bottom  row  is  expanded  into  the  following  equation; 

hdfloXo  "FEbO^B  "^^BB^B 

From  Eq.  (28),  the  base  acceleration  is  obtained  as  follows; 

Xg  ~^Bb1^BO^O  "F^BB^bJ  (^^) 

From  Eq.  (2S),  the  top  row  is  expanded  to  obtain: 

^00^0  "Fb^oB^B  "Fl^oo^o  "FI^oB^B 

Substituting  Eq.  (94)  into  Eq.  (74)  and  simplifying,  the  equation  of  motion  for  the  internal 
coorchnates  as  a  function  of  prescribed  base  displacement  is  obtained  as  follows; 

[b^oo  “bdQgMggMpQ Jcq  +|Kqq  "MggMggKooJxQ  =Fq  ~{Kob  "FMoBbdBBKggJtg  (95) 
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N(Mwe  that  the  matrix  that  pre-muhiplies  the  acceleration  term  has  units  of  mass,  and  the 
matrix  that  pre-multiplies  the  displacement  term  has  units  of  stiflfiiess.  The  matrix  that 
pre-multiplies  the  base  displacem<mt  term  has  units  of  stiffiiess.  Therefore,  all  terms  of  this 
equation  of  motion  are  dimensionally  consstent  with  units  of  force. 

The  same  nrodal  decomposition  pro<^ures  described  in  the  base  acceleration 
formulation  apply  to  the  base  displacement  formulation.  Since  the  di^lacemoits  at  the 
base  are  prescribed,  the  base  coordinates  are  no  longer  degrees  of  freedom.  This  means 
that  there  are  no  rigid  body  modes  associated  with  the  base  displacement  formulation.  As 
a  result,  Eq.  (89)  is  used  to  obtain  the  response  of  all  coordinates.  As  in  the  base 
acceleration  formulation,  the  time  hi^ory  of  the  base  di^lacement  is  taken  to  be  simple 
harmonic,  and  is  represented  by; 

Numerical  convergence  examples  are  provided  in  the  following  section.  Both  examples  are 
similar  to  the  examples  presented  in  the  prescribed  base  acceleration  problem. 

1.  Tip  Deflection  Calculation 

The  mast  and  antenna  system  was  subjected  to  a  base  excitation  where  the  time 
history  of  the  displacement  of  the  base  coordinates  was  prescribed.  The  excitation  was 
performed  at  a  frequency  which  corresponded  approximately  to  mode  S  (47.26  Hz)  of  the 
mast/antenna  system.  Since  the  system  was  modeled  without  damping,  a  frequency  which 
corresponded  exactly  to  a  natural  frequency  of  the  S3rstem  could  not  be  prescribed.  The 


excitation  fi^uency  is  equivalent  to  the  mode  S  natural  frequency  to  within  2  decimal 
places.  As  in  the  base  acceleration  problem,  the  frequency  was  arbitrarily  selected.  The  tip 
deflection  of  the  anterma  was  calculated  and  percent  error  in  tip  deflection  was  plotted 
versus  percent  of  available  component  modes  retained  (see  Figures  18  and  19).  The 
calculations  were  performed  twice.  Mast  modes  were  truncated  in  the  first  calculation, 
and  anterma  modes  were  truncated  in  the  second  calculation.  All  three  procedures  yidded 
accurate  results  in  both  the  mast  and  anterma  truncation  tests.  When  the  mast/anterma 
system  was  excited  at  mode  5,  the  Craig-Chang  and  MacNeal  procedures  provided  very 
similar  results  for  both  the  mast  truncation  and  antenna  trimcation.  As  was  determined  in 
the  tip  deflection  calculations  of  the  base  acceleration  problem,  the  Craig-Bampton 
procedure  yielded  the  best  results  using  fewer  modes  than  did  the  Craig-Chang  and 
MacNeal  procedures  in  both  truncation  tests.  Hiis  again  could  possibly  be  due  to  the 
combination  of  fixed  inter&ce  normal  modes  and  static  constraint  modes  providing  a 
better  representation  of  prescribed  base  displacement  than  the  combination  of  free 
interface  normal  modes  and  residual  fl«dbility  modes.  However,  all  three  provided 
accurate  results  and  at  a  cheaper  cost  than  the  FE  model. 
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Figure  18:  Percent  Error  in  Antenna  Tip  Deflection  Plotted  Versus  the 
Percent  of  Available  Mast  Modes  Retained  (Forcing  Frequency:  47.26  Hz) 
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Figure  19:  Percent  Error  in  Antenna  Tip  Deflection  Plotted  Versus  the 
Percent  of  Available  Antenna  Modes  Retained.  (Forcing  Frequency:  47.26  Hz) 
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2.  Moment  and  Shear  Cakulation 

Percent  error  in  shear  and  moment  at  the  mast  and  antenna  connectic  i  were 
plotted  versus  the  percent  of  available  component  modes  retained  (see  Figures  20-23). 
The  base  excitation  was  conducted  again  at  47.26  Hz  Similar  results  were  obtained  in  the 
shear  and  moment  calculations  as  were  obtained  in  the  tip  deflection  calculations.  The 
Craig-Bampton  procedure  provided  more  accurate  resuhs  using  fewer  modes  than  the 
Craig-Chang  and  MacNeal  methods.  The  convergence  rate  of  the  Craig-Chang  and 
MacNeal  methods  were  almost  identical  in  the  mast  mode  truncation.  When  the  mast  was 
excited  at  mode  S,  all  three  yielded  excelloit  results  in  the  antenna  mode  truncation 
analysis.  However,  the  Craig-Bampton  procedure  converged  more  quickly  than  the  other 
methods.  It  appears  that  the  combination  of  fixed  interfiice  normal  modes  and  static 
constraint  modes  lead  to  a  higher  rate  of  convergence  in  determf  >ing  tip  deflection  and 
antenna/mast  shear  and  moment  calculations  Although  all  three  procedures  initially  had  a 
higher  percentage  error  when  truncating  mast  modes,  than  vdien  truncating  antennae 
modes,  they  all  provide  as  accurate  if  not  more  accurate  results  than  the  antenna 
tnmcation  at  less  cost  in  terms  of  computations.  It  is  also  noteworthy  to  compare  the 
computational  cost  in  retaining  mast  modes  versus  retaining  antenna  modes  in  predicting 
accurate  system  response.  For  example,  in  the  shear  calculation  when  the  mast  was 
excited  at  mode  5,  the  Craig-Chai^  procedure  yielded  a  percentage  error  of  2.1%  while 
using  1.66  *  10^  FLOPS  dunng  the  mast  mode  truncation  test.  During  the  antenna  mode 
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tnincation  test,  the  Craig-Chang  procedure  yielded  a  percentage  error  of  2.03%  u^g 
1.89  10*  FLOPS. 

Since  the  impact  of  shock  waves  on  the  mast/antoina  system  are  typically  of  low 
frequency,  and  as  can  be  seen  from  the  results  of  the  prescribed  base  acceleration  and  base 
displacement  examples  provided,  it  is  recommended  that  the  Craig-Bampton  component 
mode  represoitation  be  used  to  synthesize  the  mast  and  amenna  system 
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Figure  21:  Percent  Error  in  Momoit  at  Mast  and  Antenna  Connection  Plotted 
Versus  the  Percent  of  Available  Mast  Modes  Retained. 

(Fordng  Frequency:  47 J6  Hz) 
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Figure  23:  Percent  Error  in  Moment  at  Mast  and  Antenna  Connection  Plotted 

Versus  the  Percent  of  Available  Antenna  Modes  Retained. 

(Forcing  Frequency:  41.7iS  Hz) 
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VL  CONCLUSIONS  AND  RECOMMENDATIONS 


The  survivability  of  shipboard  combat  systons  equipment  is  paramount  to  the  warfare 
fighting  capability  of  the  ship  and  her  crew.  Should  a  fire  control  radar,  or  any  vital 
topside  cond>at  systems  equipmoit  fiul  as  a  result  of  an  induced  shock  wave,  the  ship's  war 
fighting  ciqnicity  would  be  crippled.  However,  using  proven  structural  dynamics 
techniques,  the  design  engineer  can  design  the  mast/antenna  system  in  such  a  maimer  as  to 
minimize  risk  of  &ilure. 

This  study  has  been  conducted  to  show  tlm  progress  that  has  been  made  thus  far  in  the 
maximization  of  condiat  survivability  of  ^{^ard  mast  and  antenna  systems.  As 
demonstrated  in  this  study,  the  mast  and  antenna  can  be  treated  as  separate  substructures, 
and  using  CMS,  can  be  assonbled  as  a  mast/antenna  system,  fix>m  ^ch  dynamic 
re^nse  to  base  excitation  can  be  calculated.  Treating  each  antenna  as  a  substructure, 
allows  the  cataloging  of  the  various  antennae.  A  selected  antenna  can  be  "plugged”  into 
various  locations  along  the  mast  until  a  suitable  dynamic  response  is  obtained.  As  has 
been  dononstrated  herein,  CMS  along  with  FE  modeling  provides  rapid  and  accurate 
results  at  a  conqnitational  cost  ngnificantly  less  than  standard  FE  modeling  alone.  The 
mast  and  antenna  model  used  in  this  study  coninsted  of  only  17  elements  which 
corre^nds  to  a  total  of  51  degrees  of  fieedom.  Although  the  resuhs  that  were  obtained 
on  this  "nnall"  model  were  accurate  and  computationally  efficient,  the  same  benefits  can 


be  expected  with  laiger  models;  models  that  will  be  used  to  represent  realistic  mast  and 
antenna  systems. 

All  three  methods  yield  results  that  are  accurate  and  more  computationally  efficient 
than  standard  FE  modeling.  However,  from  the  results  obtained,  it  is  strongly 
recommended  that  the  Craig-Bampton  component  mode  representation  be  used  to 
synthesize  the  mast  and  antenna  system.  Since  the  Craig-Bampton  procedure  yidded 
more  accurate  results  while  using  fewer  component  modes  than  the  other  methods,  the 
Craig-Bampton  procedure  is  the  substructure  coupling  formulation  of  choice  due  to  good 
accuracy  and  ease  of  implementation.  It  is  strongly  suggested  that  the  following 
recommendations  be  implemented  when  computing  the  dynamic  response  of  a  shipboard 
mast  and  amenna  system; 

•  Modify  the  existing  FE  code  such  that  the  mast  and  antenna  system  can  be 
modeled  with  6  DOF/node.  This  type  of  modeling  will  allow  for  out  of  plane 
dynamic  analysis.  In  addition  to  accomodating  a  more  general  analysis  of  the  mast 
and  antenna  system,  model  the  mast  and  antennae  with  other  types  of  elements 
such  as  shell  and  plate  elements  in  addition  to  the  existing  beam  elements. 

•  Use  CMS  in  conjunction  with  FE  modeling.  In  particular,  it  is  recommended  that 
the  Craig-Bampton  procedure  be  used  as  the  substructure  coupling  method. 
Additionally,  if  another  software  package  which  is  tailored  to  FE  analysis  be  used 
to  conduct  the  dynamic  analysis,  the  NASTRAN  superelement  scheme  contains  the 


89 


Craig-Bampton  component  mode  representation  as  a  solution  path  to  dynamic 
analyses. 

In  this  study  it  has  been  demonstrated  that  the  mast  and  the  various  antennae  can  be 
modeled  separately  as  substructures  using  standard  FE  modeling.  Using  component  mode 
synthesis,  the  substructures  can  be  assembled  into  a  system  and  dynamic  response 
comiHited  accurately  at  a  computational  cost  that  is  less  than  standard  FE  modding  alone. 
Implementing  the  above  recommendations  will  lead  to  the  maximization  of  topside  combat 
systems  equipmoit  survivability. 
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APPENDIX  A 


%  LT  Lynn  James  Petersen 
%  Naval  Postgraduate  School 

%  Mast-Antenna  Survivability;  Structural  Dynamic  Design  Analysis  by  Component  Mode 
%  Synthesis 

%  Main  Program:  Progsynt.m 

%  This  program  numerically  constructs  a  FE  model  of  the  either  the“mast”,  the 
%  “antenna”  as  components  or  the  mast  and  antenna  as  a  system.  The  user  specifies  the 
%  desired  configuration 

%  To  begin  the  program,  the  user  specifies  the  desired  structure  (i.e.  the  mast,  or 
%  antenna,  or  mast/antenna  system.) 

p  =  input('enter  which  structure  is  to  be  assembled:  enter  1  if  mast/antenna,... 
enter  2  if  mast,  enter  3  if  antrana'); 
ifp  =  =  1 

datdabe  %  data  file  for  the  mast/antenna  system 
els«fp==2 

data_mst  %  data  file  for  the  mast 
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else 


(lata_aiit  %  data  file  for  the  antenna 
end; 

ilops(0)  %  reset  the  FLOPS  count 

%  A.  Calculation  of  the  number  of  elements,  nodes,  length  of  each  element,  radius  of 
gyration  and  orioitation. 

A1  =  size(con);  %  “con”  is  the  connectivity  matrix 
numel  =  Al(l);  %  number  of  elemotts 
Cl  =  size(coord);  %  “coord”  is  the  matrix  of  coordinates 
nodes  =  C(1);  %  number  of  nodes 

for  i  =  Irnumel; 

IC  =  con(i,l); 

ID  =  con(i,2); 

%  determine  the  length  of  the  respective  beam  element 

l(l,i)=sqrt((coord(ID,l)-coord(IC,l))^2+(coord(ID,2)-coord(IC,2))^2); 
DX(i)=coord(ID,  1  )-<x>ord(IC,  1 ); 

DY(i)=coord(ID,2)-coord(IC,2); 
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%  determine  the  orientation  of  the  beam  element 


if  DX(i)>=  0  &  DY(i)>=  0; 
theta(  1  ,i)=acos(DX(i)/l(  1  ,i)); 

elseif  DX(i)<0  &  DY(i)>=0; 

theta(  1  ,i)=acos(D  Y(i)/1(  1  ,i))+pi/2; 

elseif  DX(i)<0  &  DY(i)<0; 

theta(  1  ,i)=acos(abs(DX(i))/l(  1  ,i))+P>; 

else 

theta(  1  ,i)=acos(abs(DY(i))/l(l  ,i))+(3  *pi/2); 

Old; 

%  calculate  the  radius  of  gyration 

r(l,i)  =  sqrt(I(l,i)/A(l,i)); 

end; 

%  B.  Computer  truncation  error  minimization 

%  The  function  trig  will  now  be  called  to  ensure  that  “exact”  values  are  sent  to  the 
%  stifihess  and  mass  matrix  function. 

[c  1  ,s  1  ]=trig(theta,numel); 
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%  C.  Assonbly  of  stifihess  and  mass  matrices; 

%  The  function  “E13mk.m”  will  be  called  to  assemble  the  elemental  mass  and  stifihess 
matrices. 

%  The  “global”  mass  and  stifihess  matrices  will  be  assembled  simultaneously 
kbl  =  [zeros(nodes  ♦  3,nodes  *  3)]; 
mbl  =  [zeros(nodes  *  3,nodes  *  3)]; 
for  i=l:numel; 

[ke,me]=el3mk(l(i),gamma(i),I(i),E(i),r(i),cl(i),sl(i)); 

v=con(i,l); 

w=con(i^); 

kbl(3*v-2:3*v,3*v-2:3*v)  =  kbl(3*v-2:3*v,3*v-2:3*v)  +  ke(l:3,l:3); 
kbl(3*v-2:3*v,3*w-2;3*w)  =  kbl(3*v-2:3*v,3*w-2:3*w)  +  ke(l:3,4:6); 
kbl(3*w-2:3*w,3*v-2:3*v)  =  kbl{3*w-2:3*w,3*v-2:3*v)  +  ke(4:6,l  :3); 
kbl(3*w-2:3'''w,3’''w-2;3*w)  =  kbl(3*w-2:3*w,3*w-2:3*w)  +  ke(4:6,4:6); 

mbl(3*v-2:3*v,3*v-2:3*v)  =  mbl(3*v-2:3*v,3*v-2:3*v)  +  me(l:3,l:3); 
mbl(3*v-2:3*v,3*w-2;3*w)  =  nU)l(3*v-2:3*v,3*w-2:3*w)  +  me(l:3,4:6); 
mbl(3*w-2;3*w,3*v-2:3*v)  =  mbl(3*w-2:3*w,3*v-2:3*v)  +  me(4:6,l:3); 
mbl(3*w-2;3*w,3*w-2:3*w)  =  mbl(3*w-2:3*w,3*w-2;3*w)  +  me(4:6,4;6); 

end; 
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%  D.  Boundaiy  conditions 

kbl([BC],:H  ]; 
kbl(:,[BC]H  ]; 
nibl([BC],:H  ]; 
mbl(:,[BC]H  ]; 

%  E.  Internal  and  base  coordinates 

if  p  =  =  1;  %  mast  and  antenna  system 

%  partition  the  mass  and  stiifiiess  matrices  into  base  and  internal  coordinates 
[KEX,MEX,FEX]=kmbe(kbl,mbl,INTM,BSM,XBM); 
kbl=KEX; 
mbl=MEX; 

%  convert  the  force  vector  into  a  product  of  mass  *  acceleration 
%  because  of  MATLAB'S  inefficiency  in  calculating  rigid  body  eigenvalues  and 
%  eigenvectors  (natural  frequencies  and  mode  shapes)  we  will  make  the  first  three  natural 
%  frequencies  0,  and  calculate  rigid  body  modes  just  as  was  done  in  the  Craig-Chang 
%  formulation 
P  =  [123]; 

F  =  [4;length(kbl)]; 

KEX2  =  [KEX(F4^)  KEX(F,P);KEX(P^)  KEX(P,P)]; 

MEX2=  [MEX(F^)  MEX(F,P);MEX(P4?)  MEX(PJ»)]; 
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KEX3=KEX(F,F); 

KEX4=KEX(F,P); 

RBMODES=  [-inv(KEX3)*(KEX4) 
eye(length(P))]; 

%  calculate  the  natural  frequencies  and  mode  shapes 
[lamex,phiex]=fmodes(K£X2,MEX2); 
for  i=l:3; 

lamex(i,l)=0;  %  make  first  three  natural  frequencies  equal  to  zero 
end; 
lamex; 

omeg=sqit(lamex)/(2  •pi) 
phiex(:,l:3)=RBMODES; 
b2=3; 

[phrb]  =  rigid(MEX2,phiex,b2); 

%  check  orthogonality 
phiex(:,l:3)=phrb; 

tu=phi©x(:,l  :3)'*MEX2*phiex(:,  1 :3); 

modstifl^hiex(:,  1  ;3)'*KEX2*phi«c(;,  1 :3); 

flops  %  determine  FLOPS  of  FE  formulation  and  solution 

save  data3.mat  kbl  mbl  coord  con  lamex  phiex  BS INT  XB  f  FEX 
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save  data3a.inat  cl  si  1  r 


elseif  p— 2  %  mast  substructure 

[KEX,MEX,FEX]=kmbe(kbl,mbl,INTM3SM,XBM); 


kbl=KEX; 

mbl=MEX; 

%  convert  the  force  vector  into  a  product  of  mass  *  acceleration 
%  note;  this  operation  needs  to  be  performed  now  prior  to  CMS 
%  fl  =f-  (-(8.9525*2*pir2*FEX);  %  low  frequency 

%  fl=f+  (219.265*2*pi)^2*FEX;  %  high  frequency 

flops  %  determine  the  flops  count  to  assemble  the  mast 


save  datal.mat  kbl  mbl  V  C  coord  con  INT  O  fl  XBM  BSM 
save  datala.mat  cl  si  r  1 


fl=f 


flops  %  determine  the  flops  count  to  assemble  the  amenna 
save  data2.mat  kbl  mbl  O  INT  coord  con  f  V  C 
save  data2a.mat  cl  si  r  1 
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APPENDIX  B 


%  LT  Lynn  Jam^  Petersen 
%  Naval  Postgraduate  School 

%  Mast-Antenna  Survivability;  Structural  Dynanuc  Design  Analysis  by  Component 
%  Mode  Synthesis 
%  Program:  Cbfbx.m 

%  Two  substructure  synthesis  with  forcing  vector  and  base  excitation 
%  This  program  uses  the  Craig-Bampton  method  for  synthesizing  structures  toget’ 

%  The  program  loads  the  data  files  into  the  program.  Once  into  the  program,  the  program 
%  win  synthesize  the  structures  together,  and  calculate  the  coupled  natural  fi^uency  and 
%  fi^  interface  mode  sh^)es 

%  note:  this  program  will  synthes^  only  two  structures  together 
%  A.  The  substructures  wUl  now  be  loaded 
load  datal  .mat; 

kl=ld)l,ml=mbl,Vl=V,Cl=C,fl=fl; 
load  data2.mat; 

k2=ld)l,m2=mbl,V2=V,C2=C,f2=f; 
flops(0)  %  reset  the  FLOPS  coum 
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%  B.  K  matrix  partitioning  to  support  calculation  of  static  constraint  niodes 
Kwl  =kl(Vl,Vl); 

Kvcl  =kl(Vl,Cl); 

Kw2  =  k2(V2,V2); 

Kvc2  =  k2(V2,C2); 

%  C.  Partition  stiffiiess  and  mass  matrices  for  synthesis; 
kl=[kl(Vl,Vl)  kl(Vl,Cl);kl(Cl,Vl)  kl(Cl,Cl)]; 
ml=[ml(Vl.Vl)  ml(Vl,Cl);ml(Cl,Vl)  ml(Cl,Cl)]; 
k2=[k2(V2,V2)  k2(V2,C2);k2(C2,V2)  k2(C2,C2)]; 
ni2=(m2(V2,V2)  ni2(V2,C2);m2(C2,V2)  m2(C2,C2)]; 

%  D.  Calculate  static  constraint  modes: 

Cmodesl  =»  [-inv(Kwl)*Kvcl 
eye(length(Cl))]; 

Cniodes2  =  [-inv(Kw2)*Kvc2 
eye(length(C2))]; 

%  E.  Calculate  the  fixed  interface  normal  modes. 
kfixl=kl(l  :length(Vl),  1  ;length(Vl)); 
nifixl=ml(l:length(Vl),l;length(Vl)); 
kfix2=k2(l  :length(V2),  1  :length(V2)); 
mfix2=m2(  1  :length(V2),  1  :length(V2)); 

[laml  ,Nmodesl]=finodes(kfixl  ,nifixl),pause 
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%  F.  Note:  the  next  aunnumd  retains  the  desired  modes  and  i^ppends  zo-os  to  the 
%  interface  coordinates: 

y  =  inputCEnter  the  number  of  fixed  imerface  modes  that  are  desired  to  be  kept') 
[Nmodesl]  =  [Nmodesl(:,l:y);zeros(3,y)]; 

%  G.  Obtain  natural  frequencies  and  mode  shapes  for  the  second  substructure 
[lam2,Nmodes2]=finodes(kfix2,mfix2),pause 
%  H.  Note:  the  next  command  retains  the  desired  nuxies  and  ^^)ends  zeros  to  the 
%  interface  coordinates  for  the  second  structure 

s  =  inputCEnter  number  of  modes  dewed  to  be  retained  for  second  structure') 
[Nmodes2]  =  [Nmodes2(:,l:s);zeros(3,s)]; 

%  I.  Obtain  the  “reduction”  transformation  matrix  fi’om  the  fixed  interftce  normal 
%  modes  and  static  constraint  modes: 

[NC_MODESl]  =  [Nmodesl  Cmodesl]; 

[NC_MODES2]  =  [Nmodes2  Cmodes2]; 

%  J.  Obtain  the  reduced  mass  and  stififiiess  matrices: 
kl_red  =  NC_MODESl'  ♦  kl  ♦  NC_MODESl; 
ml_red  =  NC_MODESl'  *  ml  ♦  NC_MODESl; 
k2_^red  =  NC_MODES2'  *  k2  ♦  NC_MODES2; 
m2_red  =  NC_MODES2'  *  m2*  NC_MODES2; 
a  ==  length(Cl); 

b  =  length(C2); 
d  =  length(Vl); 
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e  =  lei^;tl^V2); 
sizel  =  length(kl_red); 


size2  =  length(k2_red); 

%  K.  Form  the  uiK»upled  stiffiiess  and  mass  matrices: 

KjuncpI  =  [kl_red  zeros(sizel,size2);zeros(size2,sizel)  k2_red]; 

Mjuncpl  =  [ml_red  zeros(si2el,size2);zeros(size2,sizel)  m2_red]; 

%  L.  Place  the  force  vector  in  the  correct  form  and  pre  muhiply  by  the  transformation 
matrix: 

fl=[fl(Vl,l);fl(Cl,l)]; 

£2=(f2(V2,l);e(C2,l)]; 

FIV  =  Nmodesl(l:d.l:y)'*fl(l:d,l); 

FlC  *  Cmodesl(l:d,l:a)'*fl(l;d,l)  +  fl(d+l:length(fl),l); 

F2V  =  Nmodes2(l:e,l:s)’*f2(l:e,l); 

F2C  =  Cmodes2(l:e,l:b)'*G(l:e,l)  +  G(e+l:length(G),l); 

FCB  =  [  FIV; 

F2V; 

F1C+F2C]; 

%  M.  Form  the  coupled  stififiiess  and  ma»  matrices: 

Kcpl  =  2eros(y+s+a,y+s+a); 
g=length(K_cpI); 
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%  structure  1;  diagonal 

K_cpl(l  :y,  1  ;y)=K_cpl(  1  :y,  1  :y)+K_uiicpl(  1  :y.  1  :y); 

%  structure  2;  diagonal 

K_cpl(y+ 1  :y+s,y+'  1  :y+s)=K_cpl(y+ 1 :3H-s,y+ 1  :y+s)+. . . 

K_uncpl(y+a+ 1  ;y+a+s,y+a+ 1  ;y+a+8); 

%  structure  1&2:  connection 

K_cpl(y+s+ 1  :y+s+a,y+s+ 1  :y+s+a)=K_cpl(y+s+ 1  :y+s+a,y+s+ 1  :y+s+a)+. . . 
K_uncpl(y+1  ;y+a,y+l  :y+a)<-K_uncpl(y+a+s+l  :y+«+a+b,y+a+s+l  :y+s+a+b); 

%  mass  matrix 
M_cpl  =  2eros(y+s+a,y+s+a); 
li=lcngth(M_cpl); 
sMk 

%  structure  1:  diagonal 

M_cpl(l  :y,  1  ;y>=M_cpl(  l  :y,  I  ;y)+M_uiiq)l(l  :y,  1  :y); 

%  structure  2:  diagonal 

M_cpl(y+1  :y+s,y+l  :y+s)=M_cpl(y+l  :y+s,y+l  :y+s)+... 

M_uncpl(y+a+l  :y+a+s,y+a+l  :y+a+s); 

%  structure  1&2:  diagonal 

M_cpl(y+s+ 1  :y+s+a,y+s+ 1  :y+s+a)=M_cpl(y+s+ 1  :y+s+a,y+s+ 1  :y+s+a)+. . . 
M_unqrf(y+1  :y+a,y+l  :y^ti)+Mjuncpl(y+a->^l  :y+a^^■H),y^ll■rt+l  :yH+a+b); 

%  structure  1 ;  off  diagonal  partitions 

Mjcpl^+s+l  ;y4-s+a,l  :y)=M_cpl(y+s+l  ;y+s+a,l  :y)+Mjuncpl(y+l  ;y+a,l  ;y); 
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M_q)l(  1  :y,y^^  1  :y+s+a)=M_q)I(  1  :y,y+s+ 1  :y+s+a)+M_unq)l(  1  :y,y+ 1  :y+a); 

%  structure  2;  off  diagonal  partitions 

M_cpl(y+s+ 1  :y+s+b,y+ 1  :y+s)=M_cpl(y+s+ 1  :y+s+b,y+ 1  :y+s)+. . . 
M_uncpl(y+a+s+ 1  :y+a+s+b,y+a+ 1  :y+a+s); 

M_cpl(y+1  :y+s,y+s+l  :y+s+b)=M_cpl(y+l  ;y+s,y+s+l  :y+s+b)+. . . 
M_uncpl(y+a+ 1  :y+a+s,y+a+s+ 1  ly+a+s+b); 

%  N.  Obtain  the  C-B  component  mode  synthesized  dgenvectors  arnl  dgoivalues: 

=  finodes(K_cpl>I_cpI); 

omeg_cb=sqrt(lam_cpiy(2*pi);  %  convert  from  (rad/sec)^2  to  Hz 
%  P.  Perform  transformation  to  support  base  notation  calculation,  clean  up  first  three 
%  natural  fiequencies,  and  first  three  modes 
P=Pei]gth(Kjcpl)-2:Ieng;th(Kjcpl)]; 

F=[l;length(K_cpI)-3]; 

KEX2=[K_cpI(F40  K_cpl(F4»);  K_cpI(P4?)  K_cpl(P,P)]; 

MEX2=[M_cpl(F,F)  M_cpI(FJ*);  M_cpl(P,F)  M_cpl(P,P)]; 

KEX3=K_cpl(F,F); 

KEX4=K_cpl(FJ»); 

Rbmodes=[-iiiv(KEX3)*(KEX4) 

eyeOength(P))]; 

[lainex,phiex]=finodes(KEX2,MEX2); 
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fori*l:3; 


laiiiex(4lH); 

end; 

phiex(;,l  :3)=RbnK>des; 
b2=3; 

[phrb]=figid(MEX2,phiex,b2);  %  perform  Gram-Schmidt  oithogonalization 
phiex(;,l:3)=phib; 

%  check  oithogonalization; 

tu=phiex(:,  1 :3)'*MEX2*phiex(:,  1 :3),pause 
tk=phiex(:,l  :3)"*KEX2*phiex(:,l  :3),pause 
%  redefine  variables 
K_cpl»KEX2; 

M_cpl=MEX2; 

lam_cpl=lamex; 

phi_cpb=phiex; 

flops  %  determine  the  FLOPS 

save  datalO.mat  K_cpl  M_cpl  FCB  a b  d  e  y  s  NC_MODESl  NC_MODES2 
lam_cpl  phi_cpl 
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APPENDIX  C 


%  LT  Lynn  James  Petersen 
%  Naval  Postgraduate  School 

%  Mast-Antenna  Survivability:  Structural  Dynamic  Design  Aiudysis  by  Componoit  Mode 
Synthesis 

%  Program:  Ccrfib2.m 

%  This  program  is  written  in  accordance  with  the  Craig-Chang  Residual  Flexibility 
%  method.  The  program  loads  two  data  files  containing  information  from  a  FE 
%  generated  substructure.  This  information  indudes:  a)  the  mass  and  stiffiiess  matrices  b) 
%  listing  of  internal  and  interfru^  coordinates  c)  any  forcii%  or  base  excitation  data 
%  needed  to  solve  a  dynamic  response  problem.  The  program  is  ready  to  load  the  two 
%  data  files  into  the  woiic  space 
clear 

loaddatal.mat  %  mast  data 

Kl=ld)l,  Ml=mbl,  01=0,  II  =  INT,  Rl=ll,  FI  =  fl 
load  data2.mat  %  antenna  data 

K2=kbl,  M2=mbl,  02=0, 12  =  INT,  R2=I2,  F2  =  f 
flops(0)  reset  the  flops  count  to  zero 
%  A .  Substructure  1 

%  Partition  substructure  1  into  interrud  and  interface  coordiruites 

lOS 


Kool  =K1(01,01); 

Korl  =  Kl(014^1); 

K1=[K1(01,01)  K1(01,I1);  K1(I1,01) 

M1=[M1(01,01)  Ml(01,ll);  M1(I1,01) 

%  B.  Calculate  natural  "free  interface  normal  modes"  for  structure  1 
[laml  J^Rmodesl]  =  fruxles(Kl,Ml) 
al  =length(01); 
bl  =  lengthen); 
a  ==  size(FRmodesl); 
b  =  a(2); 
for  i=l 

laml(i,l)=0;  %make  first  three  “zero”  natural  frequencies  equal  zero  to 
%  corre^nd  to  R.B. 
end; 

%  C.  Prompt  the  user  for  the  number  of  nornol  modes  dewed  to  be  kept 
%  note:  the  number  of  total  modes  equals  the  sum  of  kept  modes  and  deleted  modes 

c  -  input('enter  the  number  of  free  interface  normal  modes  desired  to  be  kept  for. . . 
structure  1  ') 

FRKmodesl  =  [Frmodesl(:,bl+l:c)];%  kept  modes 
FRDmodesl  =  [Frmodesl(:,c+l:b)];  %  deleted  modes 
lamla=diag(laml); 
lanikl=lamla(bl+l:c,bl+l;c); 
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lamdmamla(c+l:b,c+l  ;b); 
ifr^anikl=iiiv(laiiikl ); 
invlamd  1  =iiiv(lamd  1 ); 

%  D.  Obtain  the  rigid  body  modes  like  constraint  modes  were  calculated 
RBmodesl  =  [-inv(Kool)*Korl 
eye(length(Rl))]; 

%  unity  mass  normalize  and  orthogonalize  the  rigid  body  modes  Note;  since  rigid 
%  body  modes  as  calculated  by  matlab  are  not  orthogonal  wit  mass  matrix,  a  "Gram 
%  Schmidt"  procedure  needs  to  be  performed 
[phib]  =  rigid(Ml,RBmodesl,bl) 

RBmodesl  =  phib; 

%  £.  Form  the  projection  matrix  [P]  in  support  of  the  inertia  relief  solution 
Pl=[eye(al+bl) .  Ml*RBmodesl*RBmodesr]; 

GSTARl=iiiv(Kool); 

GBSTARl  =  zeros(al+bl,al+bl); 

GBSTARl(l:al,l:al)  =  (»STARl(l;al,l:al)  +  GSTARl; 

FAl  =  [zeros(al,bl);eye(bl,bl);]; 

%  F.  Obtain  the  inertia  rdief  residual  flexibility  modes 
GFLEXl  =Pr*(®STARl*Pl; 

%  as  calculated  fit>m  kept  normal  modes 
%  IRAmodesl  =  [GHFLEXI  -  FRKmodesl*invlamkl*FRKmodesr]*FAl; 

%  as  calculated  fi-om  deleted  normal  modes 


107 


%  IRAmodesl  =  (FRDmodesl  *invlam<il  •FRDino<lesr)*FAl ; 

%  RmodeslA-IRAinodesl; 

%  Note:  because  this  code  is  tailored  to  rigid  body  for  both  structures  an  inertia 
%  relief  solution  is  bdng  performed  for  both  structures  the  following  code  would 
%  apply  if  the  first  structure  was  fully  restrained 

%  Partition  the  normal  modes  into  kept  modes  and  deleted  modes 
%  FRKmodesl  =  [FRmodesl(;,l:c)]; 

%  FRDmodesl  =  [FRmodesl(:,c+l:b)]; 

%  Diagonalize  and  partition  into  kept  and  deleted  natural  frequences 
%  lamla==diag(laml); 

%  lamkl=lamla(l:c,l;c); 

%  lamdmamla(c+l:b,c+l:b); 

%  invlamkl  =inv(lamkl); 

%  invlamdl  =inv(lanidl); 

%  D.  Obtain  the  residual  flexibility  modes 
%  Obtain  the  flexibility  matrix 

%  Gl=iiiv(Kl); 

%  FAl  =  [zeros(al,bl);eye(bl,bl)]; 

%  Obt^  the  residual  flexibility  modes  from  the  kept  modes 
%  Rmodes  =  (G1  -  FRKmodesl *invlamkl  *FRKmodesr)*Fa; 

%  from  the  deleted  modes 

%  Rmodes  =  (FRDmodesl  •in^iamdl  ♦FRDmodesr)*FAl ; 
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%  E.  Calculate  the  reduced  mass  and  stiffiiess  matrices  for  structure  Itransformation 


%  matrix  (SYSMODESl):  [T] 

%  SYSmodesl  =  [FRKmodesl  Rmodes]; 

%  KSYSl  =  SYSmodesl'*Kl*SYSmodesl 

%  MSYSl  =  SYSmodesr*Ml  •SYSmodesl 

%  n.  Substructure  2 

%  A.  Partition  substructure  one  into  internal  and  inter&ce  coordinates 
Koo2  =  K2(02,02); 

Kor2  =  K2(02,R2); 

K2=[K2(02,02)  K2(02412);  K2(R2,02)  K2(R2^)]; 
M2=[M2(02,02)  M2(02412);  M2(R2,02)  M2(R2^)]; 

%  B.  Calculate  natural  "fiiee  inter&ce  normal  modes"  for  structure  I 
[lani2  Jl(niodes2]  =  finodes(K2,M2) 
a2  =  length(02); 
b2  =  length(I2) 
d  =  size(FRmodes2); 
e  =  d(2); 
for  i=l;b2; 

lam2(i,l)=0; 

end; 
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%  C.  Prompt  the  user  for  the  number  of  normal  modes  desired  to  be  kept. 

f  =  inputCenter  the  number  of  fioe  into&ce  normal  modes  desired  to  be  kept  for 
structure  2  ') 

FRKmodes2  =  [FRmodes2(:,b2+l:0]; 

FRDmodes2  =  [FRmodes2(;,fH:e)l; 

Iam2a=diag(lam2); 

Iamk2=lam2a(b2+1  :f,b2+l  :£); 

Iamd2=lam2a(f4-1  :e,f*-l  :e); 

invlamk2=inv(lamk2); 

invlamd2=inv(lamd2); 

%  O.  Obtain  the  rigid  body  modes  as  constraim  modes  would  be  obtained 
RBmodes2  =  [-inv(Koo2)*Kor2 
eye(length(R2))]; 

%  unity  mass  normalize  and  orthogonalize  the  ri^d  body  modes 
[phib]=  rigid(M2,RBmodes2,b2)  %  Gram-Schmidt 
RBmodes2  =  phrb; 

%  E.  Form  the  projection  matrix  [p]  in  support  of  the  inertia  relief  solution 
P2=[eye(a2+b2)  -  M2*RBmodes2'''RBmodes2']; 

GSTAR2  =  inv(Koo2); 

GBSTAR2  =  zCTos(a2+b2,a2+b2); 

GBSTAR2(l:a2,l:a2)  =  GBSTAR2(l;a2,l:a2)  +  GSTAR2; 

FA2  =  [zeros(a2,b2);eye(b2,b2);]; 
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%  F.  Obtain  the  inertia  i^dTresidual  flexibility  modes... 

GFLEXl  =  P2’*C®STAR2*P2; 

%  ...fixnn  the  kqrt  modes 

%  IRAmodes2  =  [(MJEX2  -  FRKmodes2*invlamk2*FRKmodes2']*FA2; 

%...or  fix>m  the  deleted  modes 

%  IRAmodes2  =  (FRDmodes2*imiamd2*FRDmodes2')*FA2; 

%  G.  Calculate  the  reduced  mass  and  stififoess  matrices  for  structure  1  transformation 
%  matrix[T] 

%  SYSmode82  =  [FRKmodesl  IRAmodes2]; 

%  KSYS2  =  SYSmodes2'*K2*SYSmodes2 

%  MSYS2=  SYSmodes2'*M2*SYSmodes2 

%  H.  Couple  the  two  substructures. 

kl  =  [IRAmodesl(al+l:al-H>l,:)  +  IRAmodes2(a2+l:a2-H>2,:)]; 
k2^nv(kl); 

mddl  =  FRDmodesl(al+l;al+bl,:)*(iiivlamdl)^2*(FRDmodesl(al+l:al+bl,:))'; 
mdd2  =  FRDmodes2(a2+l:a2-H)2,:)*(in\damd2>^2*(FRDniodes2(a2+l;a2+b2,:))'; 
ml  =k2*(mddl+mdd2)*k2; 

%  Rdniild  the  kqrt  normal  modes  to  indude  rigid  booy  modes  and  dastic  modes  and 
%  kqjt  natural  fiequendes 

FRKmodesl  =  [RBmodesl  FRKmodesl]; 

FRKmodes2  =  [RBmodes2  FRKmodes2]; 
lamkl=^amla(l  ;c,  1  :c); 
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Iaink2=lain2a(  1  ;f,  1  :f); 

%  I.  Build  the  system  mass  and  stifihiess  matrices 
MSYST  =  ZCTOs(c+f,c+f); 

MSYST(l:c,l:c)  =  MSYST(l:c,l:c)  +  eye(c,c)+.. 

(FRKmodesl(al+l  :al+bl,;)')*ml  ♦FRKmodesl(al+l  :al+bl,;); 
MSYST(l:c,c+l:c+f)  =MSYST(l;c,c+l;c+f)  +... 

(-FRKmodesl(al+l  :al+bl,:))'*ml*FRKmodes2(a2+l  :a2+b2,:); 
MSYST(c+l:c+f;i;c)  =  MSYST(l:c,c+l:c+£)'; 

MSYST(c+l  :c+f,c+ 1  :c+f)  =  MSYST(c+l  :c+f,c+l  :c+f)+eye(f,f)+. . . 

(FRKmodes2(a2+l  ;a2+b2,;)')*ml  ♦FRKmodes2(a2+l  :a2+b2,:); 

KSYST  =  zeros(c+f,c+£); 

KSYST(l:c,I:c)  =  KSYST(l:c,l;c)  +  lamkl  + 

(FRKmodesl(al+l  ;al-H)l,:)’)*k2*FRKmodesl(al+l  :al+bl,:); 
KSYST(l:c,c+l:c+f)  =  KSYST(l:c,c+l;c+f)  + 

(-FRKmodesl(al+l  :al+bl,;))'*k2*FRKmodes2(a2+l  :a2+b2,:); 
KSYST(c+l:c+f,l:c)  =  KSYST(l:c,c+l:c+f)'; 

KSYST(c+l;c+f,c+l:c+f)  =  KSYST(c+l:c+f,c+l:c+f)  +  lamk2  + 

(FRKmodes2(a2+l  ;a2+b2,:)')*k2*FRKmodes2(a2+l  :a2+b2,:); 
%  J.  DETERMINE  THE  CRAIG  CHANG  FORCE  VECTOR 
SYSmodesl  =  [FRKmodesl  IRAmodesl]; 

SYSmodes2  =  [FRKmodes2  IRAmodes2]; 

F1=[F1(01,1);F1(I1,1)]; 
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F2»[F2(02.1);F2(I2,1)]; 
FCCl  *  SYSinodesl'*Fl; 


FCC2*SYSiiiodes2'*F2; 

%  I^ermine  the  second  tnuisformation  inatrix 

%  Note:  since  the  final  system  mass  and  stiffiiess  matrices  partitions  are  known,  this 
%  transfi^rmation  matrix  need  not  be  used  to  build  the  system  mass  and  stifiBiess  matrices 
T2  =  [-kl*FRKmodesl(aH-l:al+bl,l:c)  kl*FRKmodes2(a2+l:a2+b2,l:f); 
kl*FRKmodesl(al+l:al+bl,l:c)-kl*FRKmodes2(a2+l:a2+b2,l:£); 
eye(c,c)  zcros(c,0; 

zeros(f,c)  cye(f,f)]; 

%  determine  the  system  Craig  Chang  force  ve^or  repartition  the  force  vector 
FCCl  A=FCCl(c+l  :c+bl,l); 

FCClB=FCCl(l:c,l); 

FCC2A=FCC2(fH:fH)2,l); 

FCC2B=FCC2(l:f,l); 

FCCl  =  [FCCl  A; 

FCC2A; 

FCCIB; 

FCC2B]; 

FCC  =  T2'*FCC1; 

%  note:  since  the  calculated  system  is  not  restrained,  the  following  commands  clean  up 
%  the  rigid  body  modes 
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P»C1:3]; 

F=[4:length(KSYST)]; 

KEX2*{KSYST(F,F)  KSYST(F,P);  KSYST(P,F)  KSYST(P^)]; 
MEX2=[MSYST(F,F)  MSYST(F^);  MSYST(P,F)  MSYST(P,P)]; 
KEX3=KSYST(F,F); 

KEX4=KSYST(F4*); 

Rbino<les=[-iiiv(KEX3)*(KEX4) 

cye(lefigth(P))]; 

[Iaiiiex,phiex]=fmodes(KEX2,MEX2); 

fbri=l;3; 

Ianiex(U)=0; 

end; 

pliiex(:,  1 :3)=Rbniodes; 
b2=3; 

[plirb]=rigid(MEX2,phiex,b2);  %Gram-Schniidt 
phiex(;,l:3)=phfb; 

tu^hiex(:,l  ;3)'*MEX2*phicx(:,l  :3),pause  %oithogonality 

tt=plBcx(:.l;3)'*KEX2*phiex(:,l:3),pausc 

fcF'sizcCFCC); 


FCC=  [FCC(4:fcl,l); 
FCC(1:3,1)]; 
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KSYST-KEX2; 

MSYST=MEX2; 

lain_ccrHainex; 

phi_ccif^hiex; 

om^_ccrf»sqit(lam_ccrf)/(2*pi); 

(Mii^»oin^_ccff 

fl<^  %  determiiie  total  FLOPS 

save  dataS.mat  Iain_ccrf phi_ccff MSYST  KSYST  FCC  SYSmodesl... 
SYSniodes2  T2  al  a2  bl  b2  f  c  XBM 
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APPENDIX  D 


%  LT  Lynn  James  Petersen 
%  Naval  Postgraduate  School 

%  Mast-Antenna  Survivability.  Structural  Dynamic  Design  Aiudysis  by  Component  Mode 
%  Synthesis 
%  Pir^nup;MrfH)2.m 

%  This  program  is  written  in  accordance  with  the  MacNeal  residual  flexibility  method. 

%  The  program  loads  two  data  flies  containit^  information  fit>m  an  FE  geMrated 
%  substructure  this  information  includes:  a)  the  mass  and  stiflfoess  matrices  b)  listing  of 
%  internal  and  interfoce  coordinates  c)  any  forong  or  base  exdtatimi  data  needed  to  solve 
%  a  dynamic  resp<Hise  (noUem.  The  program  is  ready  to  load  the  two  data  files  into  the 
%  work^Moe 

loaddatal.mat  %  mast  data 

Kl=kbl,  Ml=mbl,  01=0,  II  =  INT,  R1  =  11,  Fl=fl 
load  data2.mat  %  antennae  data 

K2=kbl,  M2=mbl,  02=0, 12  =  INT,  R2  =  12,  F2=f 
flops(0);  %  zero  out  the  FLOPS  count 
*/•  I.  Substructure  1 

%  A.  Partition  substructure  one  into  intonal  and  interfiice  coordinates 
Kool=Kl(01,01); 


116 


Korl»Kl(014ll); 

K1»[K1(01,01)  K1(01,I1);K1(I1,01)  K1(11,I1)]; 

MHM1(01,01)  M1(01,II);M1(I1,01)  M1(11,I1)]; 

%  B.  Cateulate  natural  "fm  interfiure  normal  modes"  for  structure  I 
[laml  JFRmodesl]  =  fniodes(Kl>fl); 
al  »length(01); 
bl » lengthen); 
a  =  size(FRmodesl); 
b»a(2); 
for  i=l;bl; 

laml(i,l)=0;  %  for  rigid  body  modes  (mly 

end; 

%  C.  Prem^rt  the  user  fin- the  number  of  normal  modes  desired  to  be  ke{>t. 

c  -  iiqwtCenter  the  number  of  free  iitterfiu^e  normal  modes  desired  to  be  kept  for... 
structure  1  ') 

%  Paititimi  die  normal  modes  into  kqn  modes  and  ddeted  modes 
FRKmodesl  =  [FRmodesl(;,bl+l:c)]; 

FRDmodesl  =  [FRniodesl(;,c+l:b)]; 
lamla^Hlia^laml); 
ianiklHamla(bl+l  :c,bl+l  :c); 

IaindlHamla(c+l  :b,c+l  :b); 
indamkl  ^nv(lamk  1 ); 
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inviaiiKl  1  =iiiv(laind  1 ); 

%  D.  Obtain  the  rigid  body  modes  as  static  constraint  modes  would  be  obtained 
RBmodesl  =  [-inv(Kool)*Korl 
eye(length(Rl))]; 

%  unity  mass  normalize  and  orthogonalize  the  rigid  body  modes 
[ldirb]  =  rigid(M14^modesl,bl);  %Gfam-Schmidt 
RBmodesl  =  phrb; 

%  E.  Form  the  projection  matrix  [P]  in  support  of  the  inertia  relief  solution 
Pl=Ceye(al+bl)  -  Ml  ♦RBmodesl  ♦RBmodesl’]; 

GSTARl  =inv(Kool); 

(ffiSTARl  =  zeros(al+bl,al-4)l); 

<»STARl(l:al,l:al)  =  GBSTARl(l:al,l:al)  +  GSTARl; 

FAl  =  [zeroa(al,bl);eye(bl,bl);); 

%  F.  Obtain  the  residual  flexibility  modes... 

GFLEXl  =Pr^GBSTARl^Pl; 

%  ...fiom  the  kept  modes 

%  IRAmodesl  =  [GFLEXl  -  FRKmodesl  ♦imiamkl  ♦FRKmodesl'J^FAl ; 

%  ...fimn  the  ddeted  modes 

%  IRAmodesl  =(FRDmodesl^invlanKil^FRDmodesr)^FAl; 

%  Note:  this  code  is  tailored  to  a  synthesis  of  a  free-fiee  structure  to  a  free-free  structure 
%  (i.e.  mast  &  antennae  to  support  base  accdoittion  formulation) 

%  the  foUowh^  code  applies  ifthe  first  structure  was  fiilly  restrained 
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%  FRKmodesl  =  [FRinodesl(:,l:c)]; 

%  FRDmodesl  =  [FRmodesl(:,c+l:b)]; 

%  Diagonalize  and  partition  into  kept  and  dieted  natural  frequencies 
%  lamla=diag(ianil); 

%  lamkl=iamla(l;c,l:c); 

%  lanKll=lamla(c+l;b,c+l:b); 

%  inviamkl  =inv(lanikl); 

%  inviamdl  =  inv(lanidl); 

%  D.  Obtain  the  residual  flexibility  modes... 

%  Obtain  the  flexibility  matrix 
%  Gl=inv(Kl); 

%  FAl  =  [zeros(al,bl);eye(bl,bl)]; 

%  ...fiom  the  kept  modes 

%  Rmodes  =  (G1  -  FRKmodesl •invlamkl*FRKmodesr)*Fa; 

%  ...from  the  deleted  modes 

%  Rmodes  =  (FRDmodesl*invlamdl*FRDmodesr)*FAl; 

%  E.  Calculate  the  reduced  mass  and  stififliess  matrices  for  stnicture  I 
%  Note:  since  the  partitions  of  the  sy^on  mass  and  stiflBiess  matrices  are  known,  the 
%  following  calculations  need  not  be  performed. 

%  Transformation  matrix  (SYSMODESl):  [T] 

%  SYSmodesl  =  [FRKmodesl  Rmodes]; 

%  KSYSl  =  SYSmodesl’*Kl  •SYSmodesl; 
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%  MSYSl  =  SYSmodesl'*Ml*SYSinodesl; 

%  II.  Substructure  2 

%  A.  Partition  substructure  two  into  internal  and  interface  coordinates 
Koo2  =  K2(02,02); 

Kor2  =  K2(02,R2); 

K2=[K2(02,02)  K2(02J12);  K2(R2,02)  K2(R2,R2)]; 

M2=[M2(02,02)  M2(02412);  M2(R2,02)  M2(R2J12)]; 

%  B.  Calculate  natural  "free  interfru:e  normal  modes"  for  structure  2 
[lam2JFRmodes2]  =  fmodes(K2>12) 
a2  =  length(02); 
b2  -  length(I2); 
d  =*  size(FRnKxles2); 
e»d(2); 
fori=l:b2 

lain2(i,l)=K);  %  rigid  body  modes  only 
end; 

%  C.  Prompt  the  user  for  the  number  of  normal  modes  desired  to  be  kept 

f » inputCenter  the  number  of  free  interface  normal  modes  desired  to  be  kept  for... 
structure  2  ') 


FRKmodes2  =  [FRmodes2(:,b2+l;f)]; 
FRDmodes2  =  [FRmodes2(;,f+l;e)]; 
Iaiii2a=diag(lam2); 
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Iaiiik2=lam2a(b2+1  ;f,b2+l  ;f); 

Iaind2=lam2a(f<'  I  ;e,f+ 1  ;e); 

iiivlaink2=Hnv(laink2); 

iiivlamd2=inv(lamd2); 

%  D.  Obtain  the  rigid  body  inodes  as  constraint  modes  would  be  obtained 
RBittodes2  =  [-inv(Koo2)*Kor2 
eye(length(R2))]; 

%  unity  mass  normalize  and  orthogonalize  the  rigid  body  modes 
[phrb]  =  rigid(M2,RBmodes2,b2);  %  Gram-Schmidt 
RBmodes2  =  phrb; 

%  E.  Form  the  projection  matrix  [P]  in  support  of  the  inertia  relief  solution 
P2=(eye(a2+b2)  -  M2*RBmodes2*RBmodes2’]; 

GSTAR2  =  iiiv(Koo2); 

GBSTAR2  =  zaos(a2+b2,a2+b2); 

(ffiSTAR2(l:a2,l:a2)  =  GBSTAR2(l:a2,l:a2)  +  GSTAR2; 

FA2  =  [zeros(a2,b2);eye(b2,b2);]; 

%  F.  Obtain  the  residual  flexibility  modes... 

aFLEX2  =  P2'*(»STAR2*P2; 

%  ...fi'om  the  kept  modes 

%  IRAniodes2  =  [GFLEX2  -  FRKmodes2*invlamk2*FRKmodes2’]*FA2; 

%  ...fiom  the  deleted  modes 

%  IRAmodes2  =  (FRDmodes2*invlamd2*FRDmodes2')*FA2; 
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%  G.  Calculate  the  reduced  mass  and  stiffiiess  matrices  for  structure  1  transformation 


%  matrix  [T] 

%  SYSmodes2  ==  [l^RKmodes2  IRAmodes2]; 

%  KSYS2  =  SYSmodes2'*K2*SYSmodes2;  %  note:  since  final  matrix  is  known, 

%  these  calculations  are  not  required 

%  MSYS2=  SYSmodes2'*M2*SYSmodes2; 

%  H.  Couple  the  two  systems. 

kl  =  [IRAmodesl(al+l:al+bl,;)  +  IRAmodes2(a2+l:a2+b2,:)]; 
k2=inv(kl); 

%  R^uild  the  kept  normal  modes  to  include  rigid  body  modes  and  elastic  modes  and  kept 
%  natural  firequencies 

FRKmodesl  -  [RBmodesl  FRKmodesl]; 

FRKmodes2  =  [RBmodes2  FRKmodes2]; 
lanikl=lamla(l  :c,  1  :c); 

Iamk2=riam2a(  1  1  ;f); 

%  I.  Build  the  system  mass  and  stifibess  matrices 
MSYST  =  eye(c+f,c+f); 

KSYST  =  zeros(c+f,c+f); 

KSYST(l:c,l:c)  =  KSYST(l:c,l:c)  lamkl  +... 

(FRKmodesl(al+l  :al+bl,:))'*k2*FRKmodesl(al+l  :al+bl,:); 
KSYST(l:c,c+l:c+£)  =  KSYST(l:c,c+l:c+£)  +... 

(-FRKmodesl(al+l  :al+bl,:)y*k2*FRKmodes2(a2+l  :a24b2,:); 
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KSYST(c+l:c+f,l;c)  =  KSYST(I;c,c+l;c4f)'; 

KSYST(c+l:c+f,c+l:c+f)  =  KSYST(c+l;c+f,c+l:c+f)  +  laink2  +... 

(FRKmodes2(a2+l  :a2+b2,:))’*k2*FRKmodes2(a2+l  ;a2+b2,:); 

%  Build  the  transformation  matrix  for  both  »ibstructures 
SYSmodesl  =  [FRKmodesl  IRAmod^l]; 

SYSmodes2  =  [FRKmodes2  IRAmodes2]; 

%  J.  Determine  the  MacNeal  force  vector 

%  partition  the  force  vector  with  internal  forces  followed  by  inter&ce  forces 
F1=[F1(01,1);F1(I1,1)]; 

F2=[F2(02,1);F2(I2,1)1; 

%  jPremultiply  the  force  vector  by  the  first  transformation  matrix 
FMl  =  SYSniodesl'*Fl; 

FM2  =  SYSmodes2’*F2; 

%  Detomine  the  second  transformation  matrix 

T2=  [-kl*FRKmodesl(al+l:al+bl,l:c)  kl*FRKmodes2(a2+l:a2+b2,l:f); 
kl*FRKmodesl(al+l:al+bl,l:c)  -kl*FRKmodes2(a2+l:a2+b2,l:f); 
eye(c,c)  2eros(c,f); 

zeros(f;c)  eyc(f,f)]; 

%  determine  the  system  MacNeal  force  vector 
%  repartition  the  force  vector 


FMlA=FMl(c+l:c+bl,l); 

FMlB=FMl(l:c,l); 
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FM2A»FM2(f<-l;fH)2,l); 

FM2B*FM2(l:f,l); 

FM1=[  FMIA; 

FM2A; 

FMIB; 

FM2B]; 

FM  =  T2'*FM1; 

%  Partition  the  system  mass  and  stiffness  matrices  to  be  suitable  for  base  excitation 
%  problem  (prescribed  acceleration) 

P==[l:3]; 

F=[4:length(KSYST)]; 

KEX2=[KSYST(F40  KSYST(F4*);  KSYST(P,F)  KSYST(PJ»)]; 
MEX2=(MSYST(FJ)  MSYST(F4»);  MSYST(P4F)  MSYST(P^)]; 
KEX3=KSYST(F40; 

KEX4=KSYST(F^); 

Rbmodes=  [-inv(KEX3)*(KEX4) 
eye(length(P))]; 

%  calculate  the  system  natural  frequencies 
[lamex,phiex]==finodes(K£X2,MEX2); 
for  i=l:3; 

iamex(i,l)=0;  %  rigid  body  modes 
end; 
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%  clean  up  the  complex  rigid  body  eigenvectors 


ptiiex(:,  1 :3)»Rbmodes; 

[phrb]=rigid(MEX2,phiex,b2);  %Gram-Schmidt 
phiex(:,l;3)=phib; 

%  check  orthogonality  and  diagonaluation 
tuF=phiex(:,  1 :3)'*MEX2*phiex(:,  1 :3),pause 
tk=phiex(:,  1  ;3)'*KEX2*phiex(:,  1  ;3),pause 

%  repartition  the  system  force  vector  to  correspond  to  mass  and  stiffiiess  matrix 
fin=si2e(FM); 
finl=fin(l,l); 

FM=[FM(4:finl,l); 

FM(1:3,1)]; 

KSYST=KEX2; 

MSYST=MEX2; 

lam_mr^lamex; 

phi_mif=phiex; 

omega_mif=sqrt(lam_mrf)/(2*pi); 

save  datad.mat  lam_mrf phi_mrf  MSYST  KSYST  FM  SYSmodesl  SYSmodes2 

T2al  a2bl  b2fcXBM 
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APPENDIX  E 


%  LT  Lynn  J.  Petersen 

%  Naval  Postgraduate  School 

%  Mast-Antenna  Survivability;  Structural  Dynamic  Design  Analysis  By  Component 

%  Mode  Synthesis 

%  Program:  Basexda.m 

%  This  program  solves  for  the  response  of  a  syston  that  is  subject  to  base  excitation. 

%  Structure  is  generated  purdy  from  FE  modeling 
load  data3. mat 

%  A.  Input  the  inimber  of  modes  dewed  to  be  used 
n  -  iiqxitCenter  the  number  of  modes  dewed:  *) 
phiex(;,n+l  :length(phiex)>=[]; 

%  B.  Place  into  modal  coordinates 
ml  =  phiex'*mbl*phiex; 
kl  =  phiex'*kbl*phiex; 

%  C.  Hme  and  frequency  parameters 
dt  =  0.001; 
tmax  =  0.4; 
num_stq)s  =  tmax/dt; 
t=0:dt:tmax; 
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oin^ii=sqrt(lamex); 

%  freq»8.9525 ; 

%  freq»219.265; 

om^  =  fm|*2*pi;  %(rad/sec) 

%D.  Convot  the  force  vector 
(<mii^)^2*FEX; 
fl=*[fl(4:leiigtli(fl),l); 

Fl*pliiex'*fl;  %modai  force  vector 
%E.  Solving  the  differential  equation 
fori=l:n; 

for  j*=l  :num_steps+l ; 

%  for  rigid  body  nuxles 
if  i<='3; 

c  I = 1  /omega*F  1  (i,  1  )*cos(omega*t(j)); 
c2=0; 

q(iJ)=-Fl  (i,  1  )/(omega^2)*sin(omega*t(j))+c  1  *t(j)+c2; 
else 

q(ij)  =  -Fl(i,l)*omega/(omegn(i)*(lamex(i)-... 
oniega''2))*sin(omegn(i)*t(j))+Fl(i,l)/(Iamex(i)-omega^2)*sin(omega*t(j)); 
end; 
end; 
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end; 

%  F.  Converting  back  from  modal  to  phydcal  coordinates 
x=»phiex*q; 
s=’size(x) 

xHx(s(I)-2:s(l),:); 

x(l:s(l)-3,:)]; 

x*=xl; 

save  data6.mat  x 


128 


APPENDIX  F 


%  LT  Lynn  JaiMs  Petersen 
%  Naval  Postgraduate  School 

%  Mast-Antenna  Survivability;  Structural  Dynamic  Design  Analysis  by  Component  Mode 
%  Synthesis 
%  Program  Basexcb.m 

%  This  program  solves  for  the  response  of  a  system  that  is  subject  to  base  excitation 
%  using  the  Craig  -Banqiton  formulation 
load  datalO.mat 

%  A.  Input  the  number  of  modes  desired 

n  =  inputCenter  the  number  of  modes  desired; ') 
phi_cpl(;,n+ 1  ;lef%th(phi_cpl)H]; 

%  B.  Place  into  modal  coordinates 
ml  =  phi_cpr*M_cpl*phi_cpl; 
kl  =  phi__cpl'*K_cpl*phi_cpl; 

Fl=phi_cpl’*FCB; 


%  C.  Unie  and  fiequency  parameters 
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dt»  0.001; 


tinax»0.4; 

mim_stq>s  =  tmax/dt; 
t »  0:dt:tinax; 
om^  « sqit(lam_cpl); 

%  fieq»8.9525; 

%  freq»219.265; 

OTi^  =  freq*2*pi;  %(rad/sec) 
i  D.  Solving  the  difierential  equation 
fi)ri®l:n; 

fi)r  j  =  l:nuni_steps+l; 

%  for  rigid  body  inodes 

if  i<*3;  %3  rigid  body  modes 

cl-l/oni^*Fl(i,l)*cos(oin^*t(j)); 

c2=0; 

q(ij)=-F  1  (i,  1  )/(oin^^2)*sin(oinega*t(j))+c  1  *t(j)+c2; 
else 

q(ij)  =  F10,l)*oin^a/(oin^o(i)*(lam_q>l(i)onKga^))*sin(om^nCi)*t(j)) 
+  Fl(i,l)/(lam_q)l(i,l)-oin^^2)*sin(oinega*t(j)); 
end; 
end; 
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end; 


%  E.  Converting  from  modal  (phys,inodal)  to  Craig-Bampton  coordinates 
xT>W_cpl*q; 

%  F.  Cmiverting  back  from  Craig-Bampton  coordinates  to  physical 
x3»(x(l:y,:);x(yfs+l  :y+s+a,:)]; 
xll-NC_MODESl*x3; 
x4*(x(yf  1  ;y+s,:);x(y+s+l  ;y+s+a,:)]; 
x21»NC_MODES2*x4; 
x21a^x21(ef  1  ;c+b,:)pc21(  1  ;c,:)l; 
save  datall. mat  xllxlla 
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APPENDIX  G 


%  LT  Lynn  James  Petersen 
%  Naval  Postgraduate  School 

%  Mast-Antenna  Survivability:  Stnictural  Dynamic  Design  Analysis  by  Component  Mode 
%  Syntheas 
%  Program;  Basexif.m 

%  This  program  solves  for  the  re^nse  of  a  system  that  is  subjected  to  base  excitation 
%  using  the  Craig-Chang  and  MacNeal  fonmilations 

%  the  following  input  prompts  the  user  as  to  vdiat  method  will  be  used  in  the  base 
%  exdtation  problem 

typeFii^Hit('enter  type  of  readual  flexibility  method,  1  if  Craig-Chang,  2  if  Macneal*) 
if  type==l 
load  dataS.mat 
else 

load  datab.mat 
end; 

%  A  Input  the  number  of  modes  desired 

n  » input('enter  the  number  of  modes  deared:  *) 
if  type==l 
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phi_ccrf(;,n+l  ;length(plii_ccrf))={]; 
else 

phi_mrf(:,n+l;length(phi_inff))  =(]; 
end; 

%  B.  Place  into  modal  coordinates 
if  typeF'=l 

ml  =  phi_ccif  *MSYST*phi_ccif, 
kl  =  phi_ccrf  *KSYST*phi_ccrf, 

FI  =  phi_ccrf*FCC; 

else 

ml  -  phi_mrf ♦MSYST*phi_mrf, 
kl  =phi_mrf*KSYST*phi_mrf; 

FI  =  phi_mrf  *FM; 
end; 

%C.  Time  and  fiequency  parameters 
dt  =  0.001; 
tmax  =  0.4; 
niim_stq)s  =  tmax/dt; 
t  =  0;dt:tmax; 
iftype“=l 

omegn  =  sqit(lam_ccrf); 
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dse 

omegn  =  sqrt(laffl_mrO; 
end; 

%  fieq=8.9525; 

%  freq=2l9.265; 

omega  =  freq*2*pi;  %(rad/sec) 

%  D.  Solving  the  differential  equation 

if  type®  =1 
for  i=l:n; 

forj  =  l:num_steps+l; 

%  for  rigid  body  modes 

if  i<»3; 

cl®l/omega*Fl(i,l)*cos(om^*t(j)); 

c2=0; 

q(y  )=-F  1  (i,  1  )/(omega^2)*sin(omega*t0))+c  1  *t(j)+c2; 
dse 

q(ij)=  -Fl(i,l)*omega/(on»gn(i)*Oa®_ccif(i)- 
om^^2))*sin(oniegn(i)*t(j))  +  Fl(i,iy(lam_ccrfl[i,l)-om^''2)*sin(omega*t(j)); 
end; 

end; 

end; 
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dse; 

for  i=l;n; 

forj  =  l:nuin_steps+l; 

%  for  rigid  body  inodes 
if  i<=3; 

cl=l/omega*Fl  (i,  l)*cos(omega*t(j)); 

c2=0; 

q(i  j)=-F  1  (i,  1  V(oinega'^2)*sin(om^*t(j))+c  1  *t0)+c2; 
else 

q(ij)  =  -Fl(i,l)*oinega/(oinegn(i)*(lain_mrf(i)- 
om^^2))*sin(oinegn(i)*t(j))  +  Fl(i,l)/(lam_mrfl[i,l)-oinega'^2)*sin(oin^*t(j)); 
end; 
end; 
end; 
end; 

%  E.  Converting  from  modai(phys,modai)  to  Craig-Chang  or  MacNeal  coordinates 
if  type==l 
x=|Ai_ccrf*q; 
else 

x=phi_inrf*q; 

Old; 
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%  F.  Converting  back  from  Craig-Chang/MacNeal  coordinates  to  physical 


%  first  convert  back  to  the  individual  substructures 

sz=aze(x); 

szl=sz(l); 

x=  [x(szl-2:szl,;); 

x(l;szl-3,:)]; 
xl  =  T2*x; 

x2  =  [xl(bl+b2+l:bH-b2+c,:);xl(l:bl,;)]; 

x3  =  [xl(bl+b2+c+l:bl+b2+c+f,;);xl(bl+l:bl+b2,:)]; 

xstrl  =  SYSmodesl*x2; 

xatr2  =  SYSniodes2*x3; 

xstr2=[xstr2(a2+ 1  ;a2+b2,:);xstr2(  1  :a2,:)]; 

if  type==l 

save  dataT.mat  xstrl  xstr2 
else 

save  dataS.mat  xstrl  xstr2 
end; 
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APPENDIX  H 


%  LT  Lynn  James  Petersen 
%  Naval  Postgraduate  School 

%  Mast-Antenna  Survivability  Structural  D^uunic  Design  Analysis  by  Component  Mode 
%  Symhesis 

%  This  appendbc  contains  an  alphabetical  listing  of  the  auxiliary  functions  that  support 
%  the  main  programs  that  are  listed  in  appendices  A-G 

%  Function;  E13mk.m 

function  [ke,me]  -  el3mk(l,gamma,lJE,r,c,s) 

%  This  function  is  called  to  build  the  elemental  and  global  mass  and  sdffiiess  matrices  for 
%  a  3  DOF/node  "FE"  problem. 

%  A.  Element  stiffiiess  matrice; 

ke(l,l)=l/r*l/r*c*c  +  12*s*s; 
ke(l,2)=l/r*!/r*c*s  -  12*c*s; 
ke(l,3)=-6*l*s; 
ke(l,4)=-l/r*l/r*c*c  -  12*s*s; 
ke(l,5)=-l/r*l/r*c*s+  12*c*s; 
ke(l,6)=-6*l*s; 
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ke(2,l)*l/r*l/r*c*s-  I2*c*s; 
kc(2^)=l/r*l/r*s*s  +  I2*c*c; 
ke(2,3)=6*l*c; 

ke(2,4)=-l/r*l/r*c*s  +  12*c*s; 

kc(2,5)=-l/r*l/r*s*s  -  12*c*c; 

ke(2,6)=6*l*c; 

kc(3,l)=-6*l*s; 

ke(3^)=6*I*c; 

kc(3,3)=4*l*l; 

ke(3,4)=6*l*s; 

ke(3,5)=-6*!*c; 

ke(3,6>=2*l*l; 

kc(4,l)=-l/r*l/r*c*c  -  12*s*s; 
kc(4,2)=-l/r*l/r*c*s  +  12*c*s; 
kc(4,3)=6*l*s; 

ke(4,4)=l/r*l/r*c*c  +  12*s*s; 
ke(4,5)=l/r*l/r*c*s  -  12*c'*s; 
ke(4,6)=6*l*s; 

kc(5,l)=^l/r*l/r*c*s+  12*c*s; 
ke(5,2)=-l/r*l/r*s*s  -  12*c*c; 
ke(5,3)  =  -6*l*c; 
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ke(5,4)=l/r*l/r*c*s  -  I2*c*s; 

ke(5,5)=l/r*l/r*s*s  +  12*c*c; 

ke(5,6)=-6*l*c; 

ke(6,l)=-6*I*s; 

kc(6,2)=6*l*c; 

ke(6,3)*2*l*I; 

kc(6.4)=6*l*s; 

ke(6,5)=-6*l*c; 

kc(6,6)=4*I*l; 

ke=E*I/l'^3*ke; 

g=386.09; 


%  B.  Element  mass  matrix 

me(l,l)=140*c*c+  156*s*s; 
me(l^)=-16*c*s; 
me(l,3)=-22*I*s; 
me(l,4)=70*c*c  +  54*s*s; 
me(l,5)=16*c*s; 
me(l,6)=13*l*s; 
me(2,l>=-16*c*s; 
me(2,2)=156*c*c+  140*s*s; 


iiie(2,3>»22*i*c; 

ine(2,4)=16*c*s; 

iiie(2,5)=54*c*c  +  70*s*s; 

iiie(2,6)=-13*l*c; 

inc(3,l)=-22*I*s; 

me(3,2)=22*l*c; 

me(3.3)=4*I*l; 

ine(3,4)=-l3*l*s; 

ine(3,5)=13*l*c; 

iiie(3,6)=-3*l*l; 

iiie(4,l)=70*c*c  +  54*s*s; 

iiie(4,2)=16*c*s; 

iiic(4,3)=-13*!*s; 

mc(4,4)=140*c*c  +  156*s*s; 

iiic(4,5)=-16*c*s; 

mc(4,6)=22*l*s; 

iiic(5,l)=16*c*s; 

inc(5,2)=54*c*c  +  70*s*s; 

iiie(5,3>=13*l*c; 

iiie(5,4)=-16*c*s; 

iiie(5,5)=156*c*c+  140*s*s; 


iiie(5,6)=-22*l*c; 

nie(6,l)=13*l*s; 

me(6,2)=-13*l*c; 

me(6,3)=-3*l*l; 

iiie(6,4)=22*l*s; 

iiie(6,5)=-22*l*c; 

nie(6,6)=4*l*l; 

me  =  gaiiiina*l/(420)*me; 

%  Function;  Fmanor 

%  fiinction[phi_normal,orth]==finanor(phi,n]ass) 

%  this  function  mass  nonnalizes  a  modal  matrix. 
a  =  size(phi); 
nummodes«a(l,2); 
phi_normal=zeros(a); 

% 

for  ii=l  inummodes; 

nK)dalmass(u)==phi(:,ii)'*mass*plii(:,ii); 
if  modalmass(ii)~=0, 

plii_norTnal(;,ii)=(l/sqrt(modalma5s(ii)))*phi(;,ii); 

else 
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phi_iioniial(:,u)=phi(;,u); 

end; 

end; 

%  do  the  ortho  caic: 

% 

orth=phi_nonnal'*inass*phi_nonnaI*  1 00; 

%  Function:  Fniodes.ni 

%  this  function  returns  a  vector  containing  nrode  freqs  (rad/sec)^2  and  a  matrix 
%  containing  the  mass 

%  normalized  mode  shapes.  The  mode  information  is  sorted  by  frequency  in  ascending 
%  order. 

%  [oni^a,phi]»finodes(k,m) 
function  [om^a,phi]=fniodes(k,m) 
a^ength(m); 

[v,d]=eig(m\k); 

[omga,index]=sort(diag(d)); 

(Mnega==’zeros(a,a); 

fori=l:a; 

om^a(i,i)=omga(i); 

end; 
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fori*l:a; 

phiteitip(:,i)=v(:,indcx(i)); 

end; 

(Mn^iPKliag(ofn^); 

[phi,orthp»finanoF(phitenip,m); 

%  Function:  Kmbe.m 

function  [KEXAfEXJ^X]=knibe(kbl,inbl,INT3S,XBM) 

%  A.  Partitioning  of  the  nuss  and  sdfihess  matrices  according  to  internal  and  base 
coordinates 

Mn*nibl(INT,INT); 

NUB  =  nibl(INT3S); 

MBI  =  nibl(BS,INT); 

MBB  »  nibl(BS3S); 

Kn  =  kbl(INT,INT); 

KlB»kbl(INT3S); 

KBl  =  lcbl(BS,INT); 

KBB  »  kbl(BS3S); 
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%  B.  Form  the  ex-mass  and  ex-stifi&iess  matrices 
MEX  >  (Mil  -  KIB*inv(KBB)*MBI); 

KEX  » (KII  -  KIB*iiiv(KBB)*KBI); 

%  C.  F<mn  the  ex-force  vector 

FEX  =  (MIB-KIB*ifiv(KBB)*MBB)*XBM; 

%  Function:  R^idm 

function  [phit)]  =  rigid(mbl,phi,b2) 

%  This  function  obtains  orthogonality  of  3  linearly  independent  but  not  orthogonal 
%  vectors  using  the  Gram-Schmidt  procedure 
if  b2=3;  %  3  rigid  body  modes 
%  1.  Obtain  “v”  vectors 

Vl=phK:,l); 

V2  =  phi(:^); 

V3  =  phK:^); 

%  2.  solve  for  alpha 

alpha  =  (V2'*mbl*Viy(Vr*rabl*Vl); 

V2T  =  V2  -  alpha*Vl; 

%  3.  solve  for  beta  and  gamma 
c(l,l)  =  V3’*mbl*Vl; 
c(2,l)  =  V3'*mbl*V2T; 
a(l,l)  =  V2r*mbl*Vl; 
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a(U)  =  Vl’*nibl*Vl; 
a(2,l)- V2r*inbl*V2T; 
a(2^)»Vr*iiibl*V2T; 

X  » iiiv(a)*c; 

beta«x(l,l); 

gniiina»x(2,l); 

V3T  =  V3  -  beta*V2T  -  gainina*Vl; 
plirb  =  [Vl  V2TV3T]; 

%  4.  unity  mass  normalize  the  phib  sht^ 
[lrfii_normal,orth]’=finanor(phrb,mbl); 
phrt) » irfu jKxmal; 
else  %2-lX)F/NODE 
\%  1.  Obtain  “v”  vectors 
VI  *plii(;,l); 

V2  =  phK:,2); 

%  2.  solve  for  alpha 

alpha  =  (V2’*mbl*Vl)/(Vr*mbl*Vl); 
V2T»V2- alpha*  VI; 

%  3.  reasign  the  rigid  body  modes 


%  4.  unity  mass  normalize  the  phrb  shapes 
[|rfii_normal,orth]=finanor(phrb^  1 ); 
phrb  ==  phi_normal; 
end; 

%  Function:  Forceda.m 

function(sh,mo)=forceda(type) 

%  this  program  calculates  the  maximum  shear  and  moment  to  a  basex  problem 

iftype==l  %  finite  element 
load  datafi.mat 
load  data3a.mat 
load  data3. mat 
datdabe 

elseif  type  =  =2%  Craig-Bampton 
load  data  ll.mat 
load  datal. mat 
lo«l  datala.mat 

dseiftype  -  =^3  %  Craig-Chang 
load  data  6.mat 
load  datal  .mat 
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load  datala-mat 


else  %MacNeal 

load  data  T.mat 
load  datal. mat 
load  dataU-mat 

%  input  location  where  shear  and  momoit  are  deared 
if  type**  1 

sm  *  iiq)ut('enter  the  element  where  the  sheo*  and  moment  are  desired') 
g*sizc(x); 
h=8(1.2); 
dseif  type*  =2 

str  *  iiq)ut(*enter  which  structure  is  desired  for  analysis) 
ifstr*“l 

load  datajnst.mat 
x*xll 
else 

load  data_ant.mat 
x*x21a 
end; 

sm  *  inputCenter  the  donent  where  the  sheer  and  moment  are  desired*) 
g*size(x); 
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h-g(U); 
dscif  type  =  =  3 

str  -  ii^)ut(‘enter  which  stnicture  is  desired  for  analysis) 
ifstr=l 

load  data_mst.mat 
x*xstrl 
else 

load  data_ant  .mat 
x=*xstr2 
end;  . 

sm  -  input('enter  the  element  where  the  she^  and  moment  are  desired') 
g=^ize(x); 
h*g(U); 


str  =  input(‘enter  which  structure  is  desired  for  analysis) 
if  str  =1 

load  datajmst.mat 
x=xstrl 
dse 

load  data_ant.mat 
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x-xstr2 


end; 

sm  ^  input('eiiter  the  element  where  the  sheer  and  moment  are  desired') 
g»si2e(x); 

[ke2,nie2]=el3nik(I(sm),gamma(sm),I(sm)»E(sm),r(sm),cl(sm),sl(sm)); 
t2=(cl(sm)  sl(sm)  0  0  0  0; 

-sl(sm)  cl(sm)  0  0  0  0; 

0  010  0  0; 

0  0  0  cl(sm)  sl(sm)  0; 

0  0  0  -sl(sm)  cl(sm)  0; 

000001]; 

%  t2=eye(6); 
if  type  =  =  1 

y3=abs(  [x(3  •cx)n(sm,  1  )-5,;)pk(3  *coiKsm,  1  )-4,;);x(3  *con(sm,  1  )- 
3,:);x(3*con(sni,2)... 

.5,:)pt(3*con(sm,2)-4,:);x(3*con(sm,2)-3,:)]); 
fiw  j=l:h; 

y4=ke2*t2*y3; 

end; 

nl  =  input  (*oiter  1  if  first  node  is  desired  or  2  if  second  node  is  desired’) 
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if  nl  ==  1 


sh“=inax(y4(2,;)) 

iTK)=inax(y4(3,;)) 

else 

sh=n»x(y4(5,:)) 

ino=nnax(y4(6,;)) 

end; 

elseif  type  <=4 
ifstr=l 

y3=abs(  [x(3*con(sin,  1  )-5,:);x(3*coi^sni,  1  )-4,:);x(3*con(sni,l)- 
3,:);x(3*con(sm,2)... 

-5,:);x(3*con(Mn,2)-4,;),x(3*con(Mn,2)-3,:)] ); 
for  j=l:h; 

y4=ke2*t2*y3; 

end; 

nl  =  input  ('enter  1  if  first  node  is  desired  or  2  if  second  node  is  desired’) 
if  nl  =  =  I 

sh=nuix(y4(2,:)) 

mo=max(y4(3,:)) 

else 

sli=max(y4(5,:)) 


ISO 


mo=inax(y4(6,:)) 


end; 

else 

y3=abs(  [x(3  •con(sin,  1  )-2,;);x(3  *con(sin,  1 )- 1  ,:);x(3  *con(sin,  1  ),;);x(3  *con(sin,2). . . 

-2,:)pc(3*coii(sin,2)-l  ,:)pt(3*con(sin,2),;)] ); 
for  j=l:h; 

y4s^c2*t2*y3; 

end; 

nl  =  input  ('enter  1  if  first  node  is  desired  or  2  if  second  node  is  desired’) 
if  nl==l 

sli=max(y4(2,:)) 

™o=niax(y4(3,:)) 

dse 

sh=raax(y4(5,;)) 

mo=niax(y4(6,:)) 

end; 

end; 

end; 
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%  Function;  Trig.m 

function  [c,s]  =  trig(theta,n) 

%  This  function  is  written  in  order  to  minimize  on  truncation  error  in  calculating  the 
%  an^e  between  beam  elements  of  finite  element  code 
forj=l:n 

if  theta(l  J)  <  0.02  &  theta(l  J)  >  6.25 
c(lj)=l; 
s(lj)=0; 

elseiftheta(lj)<  1.58  &  theta(l  j)  >  1.56 
c(lj)  =  0; 

s(ld)=l; 

elseiftheta(l  j)>  3.11  &theta(lj)<3.17 

c(lj)  =  -l; 

s(lj)  =  0; 

elseif  theta(l  J)  >  4.68  &  theta(l  J)  <  4.74 
c(lj)  =  0; 

s(lj)  =  -l 

else 

c(lj)=cos(theta(lj)); 

s(lj)=sin(theta(lj)); 

end; 

end; 
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